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The exchange bias (EB) effect was discovered 60 years ago by Meiklejohn and Bean. Meanwhile the 
EB effect has become an integral part of modern magnetism with implications for basic research 
and for numerous device applications. The EB effect was the first of its kind which relates to 
an interface effect between two different classes of materials, here between a ferromagnet and an 
antiferromagnet. Here we review fundamental aspects of the exchange bias effect. 
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I. INTRODUCTION 

The exchange bias (EB) effect was discovered 60 years 
ago by Meiklejohn and Bean (|92p . Meanwhile the EB 
effect has become an integral part of modern magnetism 
with implications for basic research and for numerous de- 
vice applications. The EB effect was the first of its kind 
which relates to an interface effect between two different 
classes of materials, here between a ferromagnet and an 
antiferromagnet. Later on the inter layer exchange cou- 
pling between ferromagncts interleaved by paramagnet 
layers was discovered (j49|, and the proximity effect be- 
tween ferromagnetic and superconducting layers was de- 
scribed (jl27l I175|l . Recent reviews on these topics can be 
found in Refs: (0 [11 [201 and also in this book. The EB 
effect manifests itself in a shift of the hysteresis loop to 
negative or positive direction with respect to the applied 
field. Its origin is related to the magnetic coupling across 
the common interface shared by a ferromagnetic (F) and 
an antiferromagnetic (AF) layer. Extensive research is 
being carried out to unveil the details of this effect, which 
has resulted in more then 600 publications in the last 
five years and since the last comprehensive reviews by 
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Nogues and SchuUer (|TUg)l . Kiwi i^, Berkowitz ([5]), and 
Stamps (fT55l) . 

An EB bilayer consists of two key elements, with rather 
different magnetic and structural properties: the ferro- 
magnetic layer and the antiferromagnetic layer. While 
the ferromagnetic layer can be studied in detail by using 
laboratory equipment like SQUID, MOKE and MFM, 
this is not the case for the magnetic interface and for the 
antiferromagnetic layer. The interface embedded in be- 
tween the F and AF layer has low volume, therefore it 
is difficult to separate its contribution from the F layer. 
Still, for the exchange bias effect the interfacial magnetic 
properties are essential for understanding the effect. For 
this purpose polarized neutron scattering and soft-xray 
magnetic scattering techniques can reveal some of the key 
magnetic properties of the interface. The AF layer has 
in principle no macroscopic magnetization, even so the 
magnetic moment of individual atoms are rather high. 
The magnetic properties of the AF materials are tradi- 
tionally studied by neutron diffraction. In thin films, 
due to the reduced AF volume available for scattering, 
this method is rather difficult to apply. Here, soft-xray 
magnetic scattering through the linear dichroism can re- 
veal information about the magnetic properties of the AF 
layer, therefore providing useful insights into the origin 
of the EB effect. 

From the application point of view the situation ap- 
pears to be less complex. The effect is being used in 
spin valves with one pinned and one free ferromagnetic 
layer which are embedded in devices such as storage me- 
dia, readout sensors, and magnetic random access mem- 
ory(MRAM). Nevertheless, robust, reliable and easy to 
control functional elements based on the exchange bias 
phenomenon require more understanding of the funda- 
mentals of the effect, which is further motivating research 
in this field. 

Because of recent experimental verifications of the ex- 
istence of interfacial layers by several groups (|54l l78| 11181 
[THIIiniliSQi 131; 136; 137; 145; 170: 172D, earlier mod- 
els of EB need to be revisited and eventually modified 
to take into account the effects, which are introduced by 
this layer. 

Here we first review some of the basic models for ex- 
change bias. We focus on numerical calculations and an- 
alytical treatment of those models which are based on the 
Stoner-Wohlfarth model (P^IMIl) . This has the advan- 
tage that analytical expressions can be derived and a nu- 
merical analysis can be much more efficiently performed 
as compared to micromagnetic simulations. It has, how- 
ever, the disadvantage that only coherent magnetization 
reversal processes are described within this formalism. 
Nevertheless, for a large fraction of experimental situ- 
ations the Stoner-Wolhfarth approach is adequate. The 
behavior of the F and AF spins during the magnetization 
reversal, as well as the dependence of the critical fields 
on the parameters of the F and AF layer are analyzed 
in detail. The Meiklejohn and Bean dHIl IHIISSI) model 
and the Mauri models are revisited and numerical 



and analytical expressions are compared. We continue 
with describing the main features of the Random Field 
model of Malozemoff (gl [HSl [Ml) and Domain State(DS) 
model (0 El [SHI [Ml nSl [m [im HHl). Then,we re- 
view the Kim and Stamp approach (!69| which focuses 
on spring-like behavior of the AF layers and coerciv- 
ity enhancement. We continue with one of the most 
recent models for exchange bias, the Spin Glass (SG) 
model p33p . Assuming a realistic state of the interface 
between the F and AF layers, the SG model describes 
well most of the important features of EB heterostruc- 
tures, including azimuthal dependence of exchange bias 
field and coercivity, AF and F thickness dependence, the 
inverse linear dependence on the lateral extension, and 
training effects. Finally we will discuss recent experi- 
ments in the light of the presented models. However, 
the main emphasis of this review is to describe the ba- 
sic models in a systematic fashion and to compare them 
with recent experimental results. 

II. STONER-WOHLFARTH MODEL 

The term anisotropy refers to the orientation of the 
magnetic moments with respect to given geometrical di- 
rections. In bulk materials the crystal axes are the ref- 
erence directions, while in thin films other reference sys- 
tems become important. In order to account for the ori- 
entation of the magnetic moments in magnetic materials, 
the minimum energy state is provided by analysis of the 
different contributing terms to the total magnetic energy: 
Zeeman term, anisotropy terms, and exchange coupling 
terms. This evaluation is performed by minimization of 
the total magnetic energy with respect to various param- 
eters. 

FcrroraagBctir flLm 
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FIG. 1 Definition of angles and vectors used in Stoner- 
Wohlfarth type model calculations. The reference direction 
of the film is along the unidirectional anisotropy. 

In the following we use the simplest possible expres- 
sion for the total magnetic energy for a ferromagnetic 
thin film and calculate the magnetic hysteresis loops. We 



3 



assume that all spins are confined in the film plane and 
that the film has a uniaxial anisotropy. The response of 
the magnetization to an applied magnetic field is then 
uniform. Therefore the spins will coherently rotate dur- 
ing the variation of the external field. The direction of 
the magnetization can be described by only one param- 
eter, namely the (9 — (3) angle defining the direction of 
the magnetization with respect to the applied field (see 
Fig. [ij. Many complexities of the magnetization reversal 
are neglected in this approach. Nevertheless, the Stoner- 
Wohlfarth (SW) model ^ [TMl [Ml, named after the 
investigators who developed it for treating the magne- 
tization reversal of a small single domain, is used with 
considerable success for various magnetic thin films and 
heterostructures . 
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FIG. 3 a) The azimuthal dependence of the normalized coer- 
cive field of a ferromagnetic film with uniaxial anisotropy. The 
curve is calculated with Eq. [s] b) The azimuthal dependence 
of the normalized irreversible switching field of a ferromag- 
netic film with uniaxial anisotropy. The curve is calculated 
using Eq. [7| 
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FIG. 2 The longitudinal (left column) and transverse (right 
column) components of the magnetization for a film with in- 
plane anisotropy. The curves are generated by numerical eval- 
uation of Eq. [3] 

The total magnetic energy per unit volume of a ferro- 
magnetic film with in-plane uniaxial anisotropy reads: 

Ev{l3) = -fioH Mpcosie - (3) + Kf sin2(/3), (1) 

where the first term is the Zeeman energy contribu- 
tion and the second term is the magnetic crystalline 
anisotropy (MCA), here assumed to have an uniaxial 
symmetry. The other parameters are: H for the ap- 
plied field, Mp for the saturation magnetization of the 
ferromagnet, Kp for the volume anisotropy constant of 
the ferromagnet, and 9 for the orientation of the applied 
magnetic field with respect to the uniaxial anisotropy di- 
rection, and P the orientation of the magnetization vector 
during the magnetization reversal. 



The minimization of the total magnetic energy with 
respect to the angle f3 and the stability equation: 



dEviP) 



= 0, 



d(3 5/32 
leads to the following equations: 



> 



— /iQ H Alp sin( 
fioH Mp cos(/3 



(3)+Kp sin(2/3) = 
) + 2Kp cos(2 /3) > 



(2) 



(3) 
(4) 



By solving Eq.[3]with the condition imposed by the Eq.|4] 
one obtains the f3 angle, which determines the longitu- 
dinal component (m|| = cos(/3 — 0)) and the transverse 
component {m± ~ sin(/3— 0)) of the magnetization. Both 
components are plotted in Fig.|2]for different in-plane ori- 
entations (6*). The evolution of the hysteresis loops for 
different angles 9 between the applied magnetic field and 
the orientation of the uniaxial anisotropy is shown in the 
left column of Fig. [2] and refiects the typical behavior of 
thin films with in-plane uniaxial anisotropy. Along the 
easy axis the hysteresis loop is square shaped and the 
transverse component is zero. When the applied field 
is perpendicular to the anisotropy axis (hard axis), the 
hysteresis loop has a linear slope, whereas the transverse 
component is ovally shaped. 

The expression for the coercive field can easily be in- 
ferred from 

Eq.|3] 



2Kf 



Mo 



Mf 



\cos0\ 



(5) 



For the hysteresis loops shown in Fig. [2] the coercive field 
follows in detail the expression above. At the position of 
the easy axis {9 ~ 0, vr) the coercive field is equal to the 
anisotropy field Ha — 2 Kp / iiq Mp , whereas along the 
hard axis {9 — ±7r/2) the coercive field is zero. Experi- 
mentally, this is an often encountered situation. For in- 
stance, polycrystalline magnetic films grown on a- plane 
sapphire substrates show such uniaxial and growth in- 
duced anisotropy due to steps at the substrate surface. 
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Aside from the coercive field dependence as a function 
of the azimuthal angle 9, another critical field can be 
recognized in Fig. |2] This is the field where the mag- 
netization changes irreversibly (i.e. where the hystere- 
sis opens). This irreversible switching field Hirr can be 
extracted by solving both Eqs. |3] and |4] Expressing 
the applied magnetic field H by its components along 
the easy and hard axis directions: H = (H^, Hy) = 
( H cos 9, H sin 9\ the solution of the system of equa- 
tions Eq. and H gives: = —Ha cos^ j3 and Hy = 
-|- Ha cos^ p. Eliminating f} from the previous two equa- 
tions we obtain the asteroid equation (^51 11531 1166p : 

\HA"' + \Hy\'"^Hll^ (6) 

Now, introducing back into the equation above the ex- 
pression for the field components we obtain for the irre- 
versible switching field the following expression (|25l 11531 
UMll : 
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FIG. 4 The asteroid curve for a film with uniaxial anisotropy. 
Two situations are depicted for finding a geometrical solution 
to the Eq[6] a) The magnetic field represented as a point PI 
lying outside the asteroid region exhibits one stable solution 
(solid line with filled arrow) and one unstable one (solid line 
with open arrow), b) A magnetic field P2 within the aster- 
oid curve exhibits four solutions(see the tangent lines) itwo of 
them are stable (solid line with filled arrow) and the other 
two are unstable(solid line with open arrow). The dotted line 
show tangents for the unstable solutions (|25l 11531 11661) . 

This field is plotted in the right panel of Fig. |3] At 
the position of the easy axis {9 — 0, tt) the irreversible 
field is equal to the anisotropy field Ha, whereas at 
9 = 7r/4, 37r/4, the irreversible field is equal to half of the 



anisotropy field {Hirr {tt / i) = Ha/2). The irreversible 
switching field can be experimentally extracted from the 
transverse components of the magnetization, whereas the 
coercive field is extracted from the longitudinal compo- 
nent of the magnetization (sec Fig. [2]). 

Fig. |4] shows the so called asteroid curve which defines 
stability criteria for the magnetization reversal(Eq. |6|. 
The asteroid method refers to an elegant geometrical so- 
lution of Eq. [l] introduced by Slonczewski (|153p . An ex- 
tended analysis can be found in Ref . (166^ . The field mea- 
sured in units of Hq appears as a point in Fig. |4j Given 
a field PI outside the asteroid curve, two solutions can 
be found by drawing tangent lines to the critical curve. 
Only one is a stable solution and is given by the tangent 
closest to the easy axis, orienting the magnetization to- 
wards the field. For fields inside the asteroid curve (P2) 
four tangents leading to four solutions can drawn. Two 
solutions are stable and the other two are unstable. The 
magnetization is stable oriented along the corresponding 
tangent (fT53 | [T66 | . 



III. DISCOVERY OF THE EXCHANGE BIAS EFFECT 

The exchange bias (EB) effect, also known as unidirec- 
tional anisotropy, was discovered in 1956 by Meiklejohn 
and Bean (PTI [Ml when studying Co particles embed- 
ded in their native antiferromagnetic oxide CoO. It was 
concluded from the beginning that the displacement of 
the hysteresis loop is brought about by the existence of 
an oxide layer surrounding the Co particles. This implies 
that the magnetic interaction across their common inter- 
face is essential in establishing the effect. Being recog- 
nized as an interfacial effect, the studies of the EB effect 
have been performed mainly on thin films consisting of a 
ferromagnetic layer in contact with an antiferromagnetic 
one. Recently, however, the lithographically prepared 
structures as well as F and AF particles are studied with 
renewed vigor. 

In Fig.|5]the original figures from Ref. ([HH [Ml) show the 
shift of a hysteresis loop of Co-CoO particles. The system 
was cooled from room temperature down to 77 K through 
the Neel temperature of CoO (T^ (CoO) = 291 K). The 
magnetization curve is shown in Fig. [5^) as a dashed line. 
It is symmetrically centered around zero value of the ap- 
plied field, which is the general behavior of ferromagnetic 
materials. When, however, the sample is cooled in a pos- 
itive magnetic field, the hysteresis loop is displaced to 
negative values (see continuous line of Fig.|5|i). Such dis- 
placement did not disappear even when extremely high 
applied fields of 70 000 Oe were used. 

In order to get more insight into this unusual effect, 
the authors studied the anisotropy behavior by using a 
self-made torque magnetometer schematically shown in 
Fig. |5ji). It consists of a spring connected to a sample 
placed in an external magnetic field. Generally, torque 
magnetometry is an accurate method for measuring the 
magnetocrystalline anisotropy (MCA) of single crystal 
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FIG. 5 a) Hysteresis loops of Co-CoO particles taken at 
77° K. The dashed line shows the loop after cooling in zero 
field. The solid line is the hysteresis loop measured after 
cooling the system in a field of 10 kOe. b) Torque curve for 
Co particles at 300 K showing uniaxial anisotropy. b) Torque 
curve of Co-CoO particles taken at 77° K showing the unusual 
unidirectioal anisotropy. d) The torque magnetometer. The 
main component is a spring which measures the torque as a 
function of the 9 angle on a sample placed in a magnetic field. 



ferromagnets. The torque on a sample is measured as a 
function of the angle between certain crystallographic 
directions and the applied magnetic field. In strong ex- 
ternal fields, when the magnetization of the sample is al- 
most parallel to the applied field (saturation), the torque 
is equal to: 

dm 

de ' 

where E{0) is the MCA energy. In the case of Co, which 
has a hexagonal structure, the torque about an axis per- 
pendicular to the c-axis follows a sin(20) function as seen 
in Fig.|5|)). The torque and the energy density can then 
be written as: 

T = -Ki sin(26') 

Ev = J Ki sm{29) dO = Ki sit? {9) + Kq, 

where Ki is the MCA anisotropy and Kq is an integra- 
tion constant. It is clearly seen from the energy expres- 
sion that along the c-axis, at 6* = and 6 = 180°, the 
particles are in a stable equilibrium. This typical case of 
a uniaxial anisotropy is seen for the Co particles at room 
temperature, where the CoO is in a paramagnetic state. 
At 77° K, after field cooling, the CoO is in an antifer- 
romagnetic state. Here, the torque curve of the Co-CoO 
system looks completely different as seen in Fig.js]:). The 
torque curve is a function of sm{9): 

T=-K.,sm{9), (8) 

hence, 

Ev^ j Ku s\n{9) d9 = -~Ku cos{9) + Kq, (9) 

The energy function shows that the particles are in equi- 
librium for one position only, namely 9 = 0. Rotating 



the sample to any angle, it tries to return to the original 
position. This direction is parallel to the field cooling 
direction and such anisotropy was named unidirectional 
anisotropy. 

Now, one can analyze whether the same unidirectional 
anisotropy observed by torque magnetometry is also re- 
sponsible for the loop shift. In Fig. [l] are shown schemat- 
ically the vectors involved in writing the energy per unit 
volume for a ferromagnetic layer with uniaxial anisotropy 
having the magnetization oriented opposite to the field. 
It reads: 

Ev = -fio H Mp cos(-/3) + Kf sin^{l3) (10) 

where H is the external field, Mp is the saturation mag- 
netization of the ferromagnet per unit volume, and Kp 
is the MCA of the F layer. The two terms entering in the 
formula above are the Zeeman interaction energy of the 
external field with the magnetization of the F layer and 
the MCA energy of the F layer, respectively. Now, writ- 
ing the stability conditions and assuming that the field 
is parallel to the easy axis, we find that the coercive field 
is: 

H, = 2Kp/tioMF (11) 

Next step is to cool the system down in an external mag- 
netic field and to introduce in Eq. [TO] the unidirectional 
anisotropy term. The expression for the energy density 
then becomes: 

Ev = -iioHMp cos(-/3) -I- Kp sin^(/3) - cos(/3) 

(12) 

We notice that the solution is identical to the previous 
case (Eq. [lO| with the substitution of an effective field: 
H' = H + Ku/Mp. This causes the hysteresis loop to 
be shifted by —K^/ iJ-o^F- Thus, Meiklejohn and Bean 
concluded that the loop displacement is equivalent to the 
explanation for the unidirectional anisotropy. 

Besides the shift of the magnetization curve and the 
unidirectional anisotropy, Meiklejohn and Bean have ob- 
served another effect when measuring the torque curves. 
Their experiments revealed an appreciable hysteresis of 
the torque (see Fig. 9 and Fig. 10 of Ref and Fig. 2 
of Ref (j9T|l ) , indicating that irreversible changes of the 
magnetic state of the sample take place when rotating the 
sample in an external magnetic field. As the system did 
not display any rotational anisotropy when the AF was 
in the paramagnetic state, this provided evidence for the 
coupling between the AF CoO shell and the F Co core. 
Such irreversible changes were suggested to take place in 
the AF layer. 

IV. IDEAL MODEL OF THE EXCHANGE BIAS: 
PHENOMENOLOGY 

The macroscopic observation of the magnetization 
curve shift due to unidirectional anisotropy of a F/AF 
bilayer can qualitatively be understood by analyzing the 
microscopic magnetic state of their common interface. 
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Phenomenologically, the onset of exchange bias is de- 
picted in Fig. [6j A ferromagnetic layer is in close contact 
to an antiferromagnetic one. Their critical temperatures 
should satisfy the condition: Tc > Tjv, where Tc is the 
Curie temperature of the ferromagnetic layer and Tjv is 
the Neel temperature of the antiferromagnetic layer. At a 
temperature which is higher than the Neel temperature 
of the AF layer and lower than the Curie temperature 
of the ferromagnet (T/y < T < Tc), the F spins align 
along the direction of the applied field, whereas the AF 
spins remain randomly oriented in a paramagnetic state 
(see Fig[6]l)). The hysteresis curve of the ferromagnet is 
centered around zero, not being affected by the proxim- 
ity of the AF layer. Next, we saturate the ferromagnet 
by applying a high enough external field Hpc and then, 
without changing the magnitude or direction of the ap- 
plied field, the temperature is decreased to a finite value 
lower then Tjv {field cooling procedure) . 
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FIG. 6 Phenomenological model of exchange bias for an AF- 
F bilayer. 1) The spin configuration at a temperature which 
is higher than Tn and smaller than Tc- The AF layer is in 
a paramagnetic state while the F layer is ordered. Its mag- 
netization curve (top-right) is centered around zero value of 
the applied field. Pannel 2): the spin configuration of the 
AF and F layer after field cooling the system through Tat of 
the AF layer in a positive applied magnetic field (Hfc)- Due 
to uncompensated spins at the AF interface, the F layer is 
coupled to the AF layer. Panel 4): the saturated state at 
negative fields. Panel 3) and 5) show the configuration of the 
spins during the remagnetization, assuming that this takes 
place through in-plane rotation of the F spins. The center 
of magnetization curve is displaced at negative values of the 
applied field by H^b- (The description is in accordance with 
Ref. ((921 [931 llMll ) 



After field cooling the system, due to exchange inter- 
action at the interface the first monolayer of the AF layer 
will align parallel (or antiparallel) to the F spins. The 



next monolayer of the antiferromagnet will align antipar- 
allel to the previous layer as to complete AF order, and 
so on (see Fig[6]-2). Note that the spins at the AF in- 
terface are uncompensated, leading to a finite net mag- 
netization of this monolayer. It is assumed that both 
the ferromagnet and the antiferromagnet are in a sin- 
gle domain state and that they will remain in this single 
domain state during the magnetization reversal process. 
When reversing the field, the F spins will try to rotate 
in-plane to the opposite direction. Being coupled to the 
AF spins, it takes a bigger force and therefore a stronger 
external field to overcome this coupling and to rotate the 
ferromagnetic spins. As a result, the first coercive field 
is higher than it used to be at T > Tjy, where the F/AF 
interaction is not yet active. On the way back from neg- 
ative saturation to positive field values (Fig|6]4)), the F 
spins require a smaller external force in order to rotate 
back (Fig[6]5)) to the original direction. A torque is act- 
ing on the F spins for all other angles, except the stable 
direction which is along the field cooling direction (uni- 
directional anisotropy). As a result, the magnetization 
curve is shifted to negative values of the applied field. 
This displacement of the center of the hysteresis loop is 
called exchange bias field, and is negative in relation to 
the orientation of the F spins after field cooling {negative 
exchange bias). It should be noted that in this simple 
description the AF spins are considered to be rigid and 
fixed to the field cooling direction during the entire re- 
versal process. 



V. THE IDEAL MEIKLEJOHN-BEAN MODEL: 
QUANTITATIVE ANALYSIS 

Based on their observation about the rotational 
anisotropy, Meiklejohn and Bean proposed a model to 
account for the magnitude of the hysteresis shift. The 
assumptions made are the following (25; 92; 109): 

• The F layer rotates rigidly, as a whole 

• Both the F and AF arc in a single domain state 

• The AF/F interface is atomically smooth. 

• The AF layer is magnetically rigid, meaning that 
the AF spins remain unchanged during the rotation 
of the F spins 

• The spins of the AF interface are fully uncompen- 
sated: the interface layer has a net magnetic mo- 
ment 

• The F and the AF layers are coupled by an ex- 
change interaction across the F/AF interface. The 
parameter assigned to this interaction is the inter- 
facial exchange coupling energy per unit area Jet 

• The AF layer has an in-plane uniaxial anisotropy 

In general, for describing the coherent rotation of the 
magnetization vector the Stoner-Wohlfarth (11591 I160p 
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model is used. Different energy terms can be added as 
needed and to best account for the quantitative and qual- 
itative behavior of the macroscopic magnetization rever- 
sal. In Fig. [7] is shown schematically the geometry of 
the vectors involved in the ideal Meiklejohn and Bean 
model. H is the applied magnetic field, which makes an 
angle with respect to the field cooling direction denoted 
by = 0, Kp and Kaf are the uniaxial anisotropy direc- 
tions of the F and the AF layer, respectively. They are 
assumed to be oriented parallel to the field cooling direc- 
tion. Mp is the magnetization orientation of the F spins 
during the magnetization reversal. It is assumed that the 
AF spins are fixed to their orientation defined during the 
field cooling procedure (rigid AF). In the analysis below 
the angle {9 = 0) for the applied field is assumed to 
be parallel to the field cooling direction. This condition 
refers to the direction along which the hysteresis loops 
is measured , whereas 9 ^ is used for torque measure- 
ments or for measuring the azimuthal dependence of the 
exchange bias field. 

IlIcjI '^IHUni^HnD ^hdri 
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FIG. 7 Schematic view of the angles and vectors used in the 
ideal Meiklejohn and Bean model. The AF layer is assumed 
to be rigid and no deviation from its initially set orientation 
is allowed. Kaf and Kp are the anisotropy of the AF layer 
and F layer, respectively, which are assumed to be parallel 
oriented to the field cooling direction. /3 is the angle between 
magnetization vector Mp of the F layer and the anisotropy 
direction of the F layer. This angle is variable during the mag- 
netization reversal. H is the external magnetic field which 
can be applied at any direction 9 with respect to the field 
coofing direction at e = (see Ref. (p2l l93l [T09)l '). 



Within this model the energy per unit area assuming 
coherent rotation of the magnetization, can be written 
as (l92] [T09l[T55] l: 

Ea = -noHMptp cos{-/3)+Kptp sin^(/3)-Jeb cos(/3) 

(13) 

where J^b [J/rn^] is the interfacial exchange energy per 
unit area, and Mp is the saturation magnetization of 



the ferromagnetic layer. The interfacial exchange energy 
can be further expressed in terms of pair exchange inter- 
actions: Eint — X^ij Jij^t'^^f^ where the summation 
includes all interactions within the range of the exchange 
couphng (43, M, ^ Hill • 

The stability condition SEa / 89 = has two types of 
solutions: one is /? = cos~^[(Je;, — ^qH Mptp)/{2Kp)] 
for H Mp tp — Jeb < 2 Kp, ; the other one is /? = 0, tt 
for iiQ H Mp tp — Jeb > '^Kp, corresponding to positive 
and negative saturation, respectively. The coercive fields 
Hci and Hc2 are extracted from the stability equation 
above for /3 — 0, tt: 



Hci 



HrO — 



2 Kptp + Jeb 
^IqMp tp 

2 Kptp — Jeb 

IJ-oMp tp 



(14) 



(15) 



Using the expressions above, the coercive field He of 
the loop and the displacement iJgfc can be calculated ac- 
cording to: 

He = and H^b = (16) 



which further gives: 



and 



He 



2Kp 



(17) 



(18) 



18 is the master formula of the EB 



The equation Eq. 
effect. It gives the expected characteristics of the hys- 
teresis loop for an ideal case, in particular the linear de- 
pendence on the interfacial energy J^b and the inverse 
dependence on the ferromagnetic layer thickness. There- 
fore this equation serves as a guideline to which experi- 
mental values are compared. In the next section we will 
discuss some predictions of the model above. 



A. The sign of the exchange bias 

Eq 



18 predicts that the sign of the exchange bias is 
negative. Almost all hysteresis loops shown in the lit- 
erature are shifted oppositely to the field cooling direc- 
tion. The positive or negative exchange coupling across 
the interface produces the same sign of the exchange 
bias field. There are, however, exceptions. Positive ex- 
change bias was observed for CoO/Co, Fea;Zni_a;F2/Co 
and Cui_a;Mna;/Co bilayers when the measuring temper- 
ature was close to the blocking temperature ([H |48l 11261 
I130j I152p . At low temperatures positive exchange bias 
was observed in Fe/FeFs (|Tn7)l and Fe/MnFa fST) bilay- 
ers. Specific of the last two systems is the low anisotropy 
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of the antiferromagnet and the antiferromagnetic type of 
coupling between the F and AF layers. It was proposed 
that, at high cooling fields, the interface layer of the an- 
tiferromagnet aligns ferromagnetically with the external 
applied field and therefore ferromagnetically with the F 
itself. As the preferred orientation between the interface 
spins of the F layer and AF layer is the antiparallel one 
(AF coupling), the EB becomes positive. Further theo- 
retical and experimental details of the positive exchange 
bias mechanism are presented in Ref . (f6T| IM) FfTjl . In the 
original Meiklejohn and Bean model the interaction of 
the cooling field with the AF spins is not taken into ac- 
count. However this interaction can be easily introduced 
in their model. The positive exchange bias could also be 
accounted for in the M&B model by simply changing the 



sign of Jeb in Eq. 13 from negative to positive. 



B. The magnitude of the EB 

Often the exchange coupling parameter Jef, is identified 
with the exchange constant of the AF layer (Jaf)- For 
various calculations a value ranging from Jaf to Jf was 
assumed. For CoO, Jaf = 21.6K = 1.86 meV dSS]). Us- 
ing this value, the expected exchange coupling constant 
Jeb of a CoO(lll)/F layer can be estimated as dHU l99)l : 



Jeb = NJaf/A = 'imJ/w?, 



where A'' = 4 is the number of Co^"*" ions at the uncom- 
pensated CoO interface per unit area A= \/(3) a? , and 

o 

a = 4.27 A is the CoO lattice parameter. With this num- 
ber we would expect for a 100 A thick Co layer, which 
shares an interface with a CoO AF layer, an exchange 
bias of: 



Heb [Oe] = „ 10^ 

Mp [kA/m] tF [A] 



(20) 



Heh — 



0.004 



1460 X 100 



10 



11 



2740 Oe 



(19) 



This exchange bias field is by far bigger than experi- 
mentally observed. So far an ideal magnitude of the EB 
field as predicted by the Eq.[T9]has not yet been observed, 
even so for some bilayers high EB fields were measured 
(see Table |l]). We encounter here two problems: first, we 
do not know how to evaluate the real coupling constant 
Jeb at the interfaces with variable degrees of complexity, 
and the second, in reality interfaces are never atomically 
smooth. The unknown interface was nicely labelled by 
Kiwi ((7?|l as "a hard nut to crack". Indeed, the features 
of the interfaces may be complex regarding the structure, 
the roughness, the magnetic properties, and domain state 
of the AF and F layers. 



TABLE I Experimental values related to Co/CoO exchange 
bias systems. The symbols used in the table are: efce-electron- 
beam evaporation, rsp-reactive sputtering, msp-magnetron 
sputtering, m6e-molecular beam epitaxy, F-ferromagnet, AF 
antiferromagnet, t^F-the thickness of the AF, tp-the thick- 
ness of the F, //eij-measured exchange bias field, //c-measured 
coercive field, Tmes-measured blocking temperature, Tmes- 
the measuring temperature, Jej-the coupling energy extracted 
from the experimental value of exchange bias field (Je6 = 
H^b (fioMptF))- 



AF 


F 


tAF 


tp 


Hcb 


H, 


Tb 


Tines 


Jeb 


Ref 






[A] 


[A] 


[Oc] 


[Oc] 


[K] 


[K] 


[mJ/m^] 




CoO (air) 


Co(rsp) 


20 


40 


-3000 


NA 




4.2 


1.75 


(97) 


CoO (air) 


Co(rsp) 


25 


27 


-2321 


3683 


180 


10 


0.91 


(ISOJ 


CoO (air) 


Co(rsp) 


25 


56 


-1073 


1751 


180 


10 


0.88 


(13^ 


CoO (air) 


Co(rsp) 


25 


87 


-675 


1315 


180 


10 


0.86 


(isni 


CoO (air) 


Co(rsp) 


25 


119 


-557 


901 


180 


10 


0.97 


(130| 


CoO (air) 


Co(rsp) 


25 


153 


-443 


789 


180 


10 


0.99 


(isni 


CoO (air) 


Co(rsp) 


25 


260 


-251 


427 


180 


10 


0.95 


(130)) 


CoO (air) 


Co(rsp) 


25 


320 


-202 


346 


180 


10 


0.94 


(13Q| 


CoO (air) 


Co(rsp) 


25 


398 


-174 


290 


180 


10 


1.00 


(1351 


CoO (air) 


Co(msp) 


33 


139 


-145 


325 




5 


0.29 


(|162|| 


CoO (air) 


Co(msp) 


33 


139 


-50 


NA 




30 


0.1 


(ITTi)) 


CoO (rsp) 


Co(rsp) 


20 


150 


-25 


295 




20 


0.055 


m 


[CoO (rsp) 


Co(rsp)]a:25 


70 


37 


-2500 


5000 




5 


13.5! 




CoO (in-situ) 


Co(ebe) 


20 


160 


-220 


330 


180 


10 


0.51 


dSI) 


CoO(lll)(mbe) 


Fe(110)(mbe) 


200 


150 


-150 


520 


291 


10 


0.4 


sec. I|114l 1134(1 



In Table [T] are listed some EB data of systems with CoO as the AF layer. We focus on experimentally de- 
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termined interfacial exchange coupling constants using 
Jeb = HebiiQ Mptp- The observed exchange couphng 
constant is usually smaller then the expected value of 
4 mJ/m? for CoO/Co bilayers by a factor ranging from 
3 to several orders of magnitude. One anomaly is seen 
for the multilayer system Co/CoO which is actually 3 
times higher then the expected value of AmJ/m?, and to 
our knowledge is the highest value observed experimen- 
tally (jl28p . Such a variation of the experimental values 
for the interfacial exchange coupling constant is motivat- 
ing further considerations of the mechanisms controlling 
the EB effect. 



for the AF spins is: a ^ 0. With this new assumption, 
the equation Eg. [Ts] reads (j9T|[T09)l : 



Ea 



-^iqH Mptp cos{e - (3) 

-Kp tp sin^(/3) -I- Kap Iaf sin^(a) 

-Jefc cos(/3-a), (22) 



where Iaf is the thickness of the antiferromagnet, and 
Kap is the MCA of the AF layer per unit area. The new 
energy term in the equation above as compared to Eq. [13] 
is the anisotropy energy of the AF layer. 



C. The 1/tF dependence of the EB field 



Eq. 18 predicts that the variation of the EB field is 
proportional to the inverse thickness of the ferromagnet: 



H. 



eh 



1 



(21) 



This dependence was subject of a large number of ex- 
perimental investigations (|109p . because it is associated 
with the interfacial nature of the exchange bias effect. 
For the CoO/Co bilayers no deviation was observed (|5T]l . 
even for very low thicknesses (2 nm) of the Co laver (|130p . 
For other systems with thin F layers of the order of sev- 
eral nanometers it was observed that the 1/tp law is not 
closely obeyed (|109p . It was suggested that the F layer 
is no longer laterally continuous p09p . Deviations from 
1/tp dependence for the other extreme when the F layer 
is very thick were observed as well (|109p . For this regime 
it is assumed that for F layers thicker than the domain 
wall thickness (500nm for permalloy), the F spins may 
vary appreciably across the film upon the magnetization 
reversal (I167p . 



D. Coercivity and exchange bias 

According to Eq. [17] the coercivity of the magnetic 
layer is the same with and without exchange bias effect. 
This contradicts experimental observations. Usually an 
increase of the coercive field is observed. 



VI. REALISTIC MEIKLEJOHN AND BEAN MODEL 



In Ref. ()93() a new degree of freedom for the AF spins 
was introduced: the AF is still rigid, but it can slightly 
rotate during the magnetization reversal as a whole as 
indicated in Fig. [8] This parameter was introduced in 
order to account for the rotational hysteresis observed 
during the torque measurements. Allowing the AF layer 
to rotate is not in contradiction to the rigid state of 
the AF layer, because it is allowed only to rotate as a 
whole. Therefore, the fourth assumption of the ideal 
M&B model in section [V] is removed. The new condition 
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FIG. 8 Schematic view of the angles and vectors used for the 
Meiklejohn and Bean model, allowing a rotation a of the AF 
layer as a whole with respect to the initially set orientation. 
Map is the sublattice magnetization of the AF layer. Kaf 
and Kp are the anisotropy of the AF layer and F layer, re- 
spectively, which are assumed to be parallel oriented to the 
field cooling direction. /3 is the angle between F magnetiza- 
tion vector Mp and the anisotropy direction of the F layer. 
This angle is variable during the magnetization reversal. H 
is the external magnetic field which can be applied at any 
direction 9 with respect to the field cooling direction at = 
(Ref. (|9ni92ll93l [Tn9|l ). 



The Eq. [22] above can be analyzed numerically by min- 
imization of the energy in respect to the a and the (3 
angles. Below we will perform a numerical analysis of 
Eq. 23 and highlight a few of the conclusions discussed 



in Ref. (j25l [911 f93|l . The minimization with respect to a 
and (3 leads to a system of two equations: 



{ 



sin(6' - (3) + sin(/3 - a) = 
sin(2 a) - sin(/3 - a) = 



where 



Jeb 



Ha Mp tp 



(23) 



(24) 
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is the value of the exchange bias field when the anisotropy 
of the AF is infinitely large, and 



R 



KaF tAF 



J. 



eb 



(25) 



is the parameter defining the ratio between the AF 
anisotropy energy and the interfacial exchange energy 
Jeb- As we will see further below, exchange bias is only 
observed, if the AF anisotropy energy is bigger than the 
exchange energy. The unknown variables a and P are 
numerically extracted as a function of the applied field 
H. Note that for clarity reasons the anisotropy of the 
ferromagnet was neglected {Kp — 0) in the system of 
equation above. As a result the coercivity, which will 
be discussed further below, is not related to the F layer 
anymore, but to the AF layer alone. Also, in order to 
simplify the discussion we consider first the case 9 = 0, 
which corresponds to measuring a hysteresis loop parallel 
to the field cooling direction. 





T 1 — J — r 



FIG. 9 Left: The phase diagram of the exchange bias field 
and the coercive fields as a function of the Meiklejohn-Bean 
parameter R. Right: Typical behavior of the antiferromag- 
netic angle a for the three difi'erent regions of the phase di- 
agram. Only region I can lead to a shift of the hysteresis 
loop. In the other two regions a coercivity is observed but no 
exchange bias field. 
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FIG. 10 Simulation of several hysteresis loops and antiferro- 
magnetic spin orientations during the magnetization reversal. 
For the simulation we used the Meiklejohn-Bean formalism. 
Top row shows three hysteresis loops calculated for different 
R ratios within region I. The graph in the last column to the 
right shows the a angle of the AF layer for the three R values. 
The middle raw shows corresponding hysteresis loops and a 
angles for R values in region II. The bottom row shows simu- 
lations for region III. Note that the scales for the a angle in 
the top panel is enlarged compared to those in the lower two 
panels. 



Numerical evaluation of the Eqs. [23] yields the angles: 

• a of the AF spins as a function of the applied field 
during the hysteresis measurement 

• /3 of the F spins which rotate coherently during 
their reversal 

The /? angle defines completely the hysteresis loop and 
at the same time the coercive fields Hd and Hc2- These 
fields, in turn, define the coercive field and the ex- 
change bias field Het (see Eq. 



16 1. The a angle infiu- 



ences the shape of the hysteresis loops when the i?-ratio 
has low values, as we will see below. For high R values 
the rotation angle of the antiferromagnet is close to zero, 
giving a maximum exchange bias field equal to . 

The properties of the EB system originate from the 
properties of the AF layer, which are accounted for by 
one parameter, the i?-ratio. We will consider the effect 
of the i?-ratio on the angels f3 and a which, as stated 
above, define the macroscopic behavior and the critical 
fields of the EB systems. 



Numerical simulations of Eq. [23' 
ratio are shown in Fig. [9] and Fig 
three physically distinct regions 



as a function of R- 
T0| We distinguish 

mi mi): 



• I. i? > 1 

In this region the coercive field is zero and the 
exchange bias field is finite, decreasing from the 
asymptotic value to the lowest finite value at 
R = 1. The AF spins rotate reversibly during 
the complete reversal of the F spins. The a angle 
has a maximum value as a function of the i?-ratio, 
ranging from approximatively zero for i? = oo to 
a = 45° at R = 1. Notice that as the maximum 
angle of the AF spins increases, a slight decrease of 
the exchange bias field is observed. When the R- 
ratio approaches the critical value of unity, the ex- 
change bias has a minimum. In the phase diagram, 
only range I can cause a shift of the magnetization 
curve. Simulated hysteresis loops for three differ- 
ent values of the i?-ratio are shown in Fig. [TO] One 
notices that not only the size of the exchange bias 
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field decreases when the i?-ratio approaches unity, 
but also the shape of the hysteresis curve is chang- 
ing. At high i?-ratios the reversal is rather sharp, 
whereas for i?-ratios close to unity it become more 
extended, almost resembling a spring-like behavior. 

• II. 0.5 < i? < 1 

Characteristic for this region is that the AF spins 
are no more reversible. They follow the F spins and 
they change direction irreversibly, causing a coer- 
cive field at the expense of the exchange bias field, 
which becomes zero. Furthermore, depending on 
the field sweeping direction, there is a hysteresis- 
like behavior of the AF spin rotation. At a critical 
angle /3 of the F spin rotation, the AF spins can- 
not withstand the torque by the coupUng to the F 
spins and they jump in a discontinuous fashion to 
another angle (jump angle). The hysteresis loops 
corresponding to this region (see Fig. 10 1 are dras- 
tically different from the previous case. The coer- 
civity shows strong dependence as function of the 
i?-ratio and they are not shifted at all. Moreover, 
the AF jump angles are clearly visible as kinks in 
the hysteresis shape during the reversal. 

• III. R < 0.5 

This region preserves the features of the previous 
one with one exception, namely that the AF spins 
follow reversibly the F spins, without any jumps. 
Therefore, no hysteresis-like behavior of the a an- 
gle is seen. The exchange bias field is zero and the 
coercive field is finite, depending on the i?-ratio. 
The hysteresis loops shown in Fig. [TO] are quite 
similar to a ferromagnet with uniaxial anisotropy. 
Within the Stoner-Wolhfarth model the resultant 
coercive field can be roughly approximated as ()25|l : 
i?c ~ 2 Kaf tAp/ (a^o Mp tp). 

It is easy to recognize that allowing the AF to rotate 
as a whole leads to an impressively rich phase diagram 
of the EB systems as a function of the parameters of the 
AF layer (and F layer). The i?-ratio can be varied across 
the whole range from zero to infinity by changing the 
thickness of the AF layer ()109p , by varying its anisotropy 
(dilution of the AF layers with non-magnetic impuri- 
ties (f60| \66\ [98|l ) , or by varying the interfacial exchange 
energy J^b (low dose ion bombardment (.23] IMl 1102(1 '). 
Recently (|17ip . an almost ideal M&B behavior has been 
observed in Ni8oFe2o/Fe5oMn5o bilayers. At high thick- 
nesses of the AF layer the hysteresis loop is shifted to 
negative values and the coercivity is almost zero, whereas 
for reduced AF thicknesses a strong increase of the coer- 
cive field is observed together with a drastic decrease of 
exchange bias. 



A. Analytical expression of the exchange bias field 

First we calculate analytically the expression of the 
exchange bias field for 6* = and Kp = 0. The exact 



analytical solution is obtained by solving the system of 
Eqs. |23]for (3 = , which leads to: 



Heb — { 



rroo 



1 







1 



R > 1 



(26) 



R < 1 



This equation retains the l/tp dependence of the ex- 
change bias field, but at the same time provides new 
features. The most important one is an additional term, 
which effectively lowers the exchange bias field when the 
R-ratio approaches the critical value of one. The R-ratio 
has three terms. One of them includes the thickness of 
the AF layer. The analytical expression for the exchange 
bias field (Eq. 26 1 predicts that there is a critical AF 
thickness Iap below which the exchange bias cannot ex- 
ist. This is 



J, 

K 



eb 
AF 



(27) 



Below this critical thickness the interfacial energy is 
transformed into coercivity. Above the critical thick- 
ness the exchange bias increases as a function of the AF 
layer thickness, reaching the asymptotic (ideal) value 
when tAF is infinite. Most recent observation of an AF 
critical thickness can be found in Rcf. ([75}. 

A similar expression as Eq. |26| was derived by Binek et 
al. ()llj) using a series expansion of Eq. 22 with respect to 
a=0 : Heb ~ H^{1 - ^) for 1/R>0. Note that close 
to the critical value of i?=l the a angle could reach high 
values up to 45°. Therefore, the series expansion with 
respect to a=0 is a good approximation for R>5. 



B. Azimuthal dependence of the exchange bias field 

In the following we consider the exchange bias field in 
region I, where it acquires non vanishing values. The 
coercive fields and the exchange bias field are extracted 
from the condition fi = + n/2 for both Hd and Hc2- 
This gives Hd = Hc2 = {—Jeb/ fJ-oMptp) cos{a{R, 9 + 
tt/2) — 0), where a{R, 9 + tt/2) is the value of the ro- 
tation angle of the AF spins at the coercive field. With 
the notation: ao = a{R, 9 + 7r/2), and using the expres- 
sion 



16 the angular dependence of the exchange bias field 



becomes: 

Heb{9) = — cos( ao - 0) 
^0 Mp tp 

The equation above can be also written as: 



Heb{e) 



Kaf tAF 
^0 Mp tp 



sin(2 ao) 



(28) 



(29) 



Interestingly, the exchange coupling parameter J^b in 
Eq.[28]is missing in Eq.|29] leaving instead an explicit de- 
pendence of the exchange bias field on the parameters of 
the antiferromagnet and the ferromagnet. The exchange 
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coupling constant and the 9 angle are accounted for by 
the AF angle a^. 



The Eq. 28 and Eq. 29 are the most general expressions 
for an exchange bias field. They include both, the influ- 
ence of the rotation of the AF layer and the influence of 
the azimuthal orientation of the applied field. Moreover, 
the anisotropy and the thickness of the AF layer are ex- 



plicitly shown in Eq. 29 To illustrate their generality we 
consider below two special cases for the Eq. [28} 

• 6* = 

In this case the hysteresis loop is measured along 
the field cooling direction {9 — 0) and Eq. 28 be- 
comes equivalent to the Eq. [26] 

• oo 

When R is very large {R ^ 1), a approximates 
zero, i.e. the rotation of the AF - layer becomes 
negligible. This is actually the original assumption 
of the Meiklejohn and Bean model. Such a condi- 
tion (-R — > oo) is approximately satisfied for large 
thicknesses of the AF layer. Then the exchange bias 
field as function of 6 can be written as (fTTl [T51 |7DJ: 



-J, 



eb 



flQ Mp t p 



cos{e) 



(30) 



In order to get more insight into the azimuthal de- 
pendence of the exchange bias field, we show in Fig. [TT] 
the normalized exchange bias field Hf,}j{9) /|i?^(0)| as a 
function of the 6 angle, according to Eq. [28] and Eq. [29] 
and for three different values of the R ratio {R— 1.1, R — 
1.5, R = 20). The ao angle (see Fig. [To|) w as obtained by 
numerically solving the system of Eqs. [23 [ For large val- 
ues of R, the azimuthal dependence of the exchange bias 
field follows closely a cos(6') unidirectional dependence. 
When, however, the R-ratio takes small values but larger 
then unity, the azimuthal behavior of the EB field devi- 
ates from the ideal unidirectional characteristic. There 
are two distinctive features: one is that at 6* = the 
exchange bias field is reduced, and the other one is that 
the maximum of the exchange bias field is shifted from 
zero towards negative azimuthal angle values. According 
to Eq. 28 this shift angle is equal to a^. In other words, 
the exchange bias field is not maximum along the field 
cooling direction. Another striking feature is that the 
shifted maximum of the exchange bias field with respect 
to the azimuthal angle 6 does not depend on thickness 
and anisotropy of the AF layer: 



MAX 



J. 



eb 



Mo Mp tp 



(31) 



Summarizing we may state that, within the Meikle- 
john and Bean model, a reduced exchange bias field is 
observed along the field cooling direction depending on 
the parameters of the AF layer {Kaf and Iaf)- How- 
ever, for i? > 1 the maximum value for the exchange bias 
field which is reached at 6* 7^ does not depend on the 
anisotropy constant {Kaf) and thickness {Iaf) of the AF 



layer. The azimuthal characteristic of the exchange bias 
allows to extract all three essential parameters defining 
the exchange bias field: J^biKAp and, Iap- 

The condition for extracting the Hd and Hc2 from the 
same (3 angle hides an important property of the magne- 
tization reversal of the ferromagnetic layer. This will be 
described next. 
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FIG. 11 Azimuthal dependence of exchange bias as a func- 
tion of the Q angle. The curves are calculated by the Eq |28| 
and Eq|29| 



C. Magnetization reversal 

A distinct feature of exchange bias phenomena is the 
magnetization reversal mechanism. In Fig. 12 is shown 
the parallel component of the magnetization TO|| = cos(/3) 
versus the perpendicular component mj^ sin(/3) for 
several R-ratios and for B — 30°. The geometrical con- 



ventions are the ones shown in Fig. [8] We see that for 
R < \ the reversal of the F spins is symmetric, similar to 
typical ferromagnets with uniaxial anisotropy. Although 
the regions 0.5 < i? < 1 and R < 0.5 exhibit different 



reversal modes of the AF spins (see Fig 10 1, there is little 
difference with respect to the F spin rotation. For both 
regions 0.5 < i? < 1 and R < 0.5 the F spins do make 
a full rotation. Similar to the uniaxial ferromagnets. At 
the steep reversal branches one would expect magnetic 
domain formation. 

When i? > 1 another reversal mechanism is observed. 
The ferromagnetic spins first rotate towards the unidi- 
rectional axis as lowering the field from positive to nega- 
tive values, and then the rotation proceeds continuously 
until the negative saturation is reached. On the return 
path, when the field is swept from negative to positive 
values, the ferromagnetic spins follow the same path to- 
wards the positive saturation. The rotation is continuous 
without any additional steps or jumps. A similar behav- 
ior was observed theoretically within the domain state 
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FIG. 12 Magnetization reversal for several values of the R 
ratio. The parallel component of the magnetization vector 
m|| = cos(/3) is plotted as a function of the perpendicular 
component of the magnetization m± — sin(/3). The reversal 
for R < 1 resembles the typical reversal of ferromagnets with 
uniaxial anisotropy. For R > 1 the reversal proceeds along 
the same path for the increasing and decreasing branch of the 
hysteresis loop. The angle 9 = 30° is chosen arbitrarily. 



model (jS] The magnetization reversal modes can be 
accessed experimentally by using the Vector-MOKE tech- 
nique (EH [In [nSl EMI) which allows to follow both the 
magnetization vector and its angle during the reversal 
process. 

D. Rotational hysteresis 

We briefly discuss again the rotational hysteresis de- 
duced from torque measurements ([Ml [55] now in 
the light of the analysis provided above. The torque mea- 
surements were carried out in a strong applied magnetic 
field H. Therefore the applied field H and the magneti- 
zation Mp can be assumed to be parallel {(3 = 0). The 
torque is given by: 

T--"^!^- Je6sin(0-aW) 

This expression differs from Eq. [8] for the ideal model 
by the rotation of the AF spins through the a angle. 
However, this does not explain the energy loss during 
the torque measurements, as observed in the experiment 
([Sjb)). The torque curve would only be a bit distorted 
but completely reversible. The integration of the energy 
curve predicts a rotational hysteresis Wrot = 0. In or- 
der to account for a finite rotational hysteresis, one can 
assume that a fraction p of particles at the F-AF in- 
terface are uniaxially coupled behaving as in region II, 
whereas the remaining fraction (1 — p) of the F-AF par- 
ticles are coupled unidirectionally, having the ideal be- 
havior as described in regime I. As seen in Fig. |9jleft), 
when the R-ratio of the uniaxial particle is in the range 
0.5 < R < 1, the AF spins will rotate irreversibly, show- 
ing hysteresis-like behavior due to a jumps indicated in 
the Fig. [9] (right). A rotational hysteresis is not expected 
for unidirectional particles with R > 1 because the AF 
structure changes reversibly with 0. With this assump- 
tion the uniaxial particles will contribute to the energy 



loss during the torque measurements, while the unidi- 
rectional particles are responsible for the unidirectional 
feature of the torque curve. This argument was used by 
Meiklejohn and Bean ^WD) when studying the exchange 
bias in core-Co/shell-CoO. A fraction p — 0.5 was in- 
ferred from the torque curves shown in Fig. [5] 

VII. NEEL'S AF DOMAIN WALL - WEAK COUPLING 

The rigid AF spin state and rigidly rotating AF spins 
concepts impose a restriction on the behavior of the 
antiferromagnetic spins, namely that the AF order is 
preserved during the magnetization reversal. Such re- 
striction implies that the interfacial exchange coupling 
is found almost entirely in the hysteresis loop either as 
a loop shift or as coercivity. Experimentally, however, 
the size of the exchange bias does not agree with the 
expected value, being several orders of magnitude lower 
then predicted. In order to cope with such loss of cou- 
pling energy, one can assume that a partial domain wall 
develops in the AF layer during the magnetization rever- 
sal. This concept was introduced by Neel (|72l I105P when 
considering the coupling between a ferromagnet and a 
low anisotropy antifcrromagnet. The AF partial domain 
wall will store an important fraction of the exchange cou- 
pling energy, lowering the shift of the hysteresis loop. 

Neel has calculated the magnetization orientation of 
each layer through a differential equation. The weak cou- 
pling is consistent with a partial AF domain wall which 
is parallel to the interface (Neel domain wall) . His model 
predicts that a minimum AF thickness is required to pro- 
duce hysteresis shift. More importantly the partial do- 
main wall concept forms the basis of further models which 
incorporate either Neel wall or Bloch wall formation as a 
way to reduce the observed magnitude of exchange bias. 



VIII. MALOZEMOFF RANDOM FIELD MODEL 

Malozemoff (1987) proposed a novel mechanism for ex- 
change anisotropy postulating a random nature of ex- 
change interactions at the F-AF interface ([5H 1551 
He assumed that the chemical roughness or alloying at 
the interface, which is present for any realistic bilayer 
system, causes lateral variations of the exchange field 
acting on the F and AF layers. The resultant ran- 
dom field causes the AF to break up into magnetic do- 
mains due to the energy minimization. By contrast with 
other theories, where the unidirectional anisotropy is 
treated either microscopically (|751 11481 I149P or macro- 
scopically (1891 [5^ [M]! , the Malozemoff approach belongs 
to models on the mesoscopic scale for surface magnetism. 

The general idea for estimating the exchange 
anisotropy is depicted in Fig. |13| where a domain wall in 
an uniaxial ferromagnet is driven by an applied in-plane 
magnetic field H ()M|) . Assuming that the interfacial en- 
ergy in one domain (cri) is different from the energy in 
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FIG. 13 Schematic side view of a F/AF bilayer with a ferro- 
magnetic wall driven by an applied field H (|84[) . 



the neighboring domain (T2, then the exchange field can 
be estimated by the equilibrium condition between the 
applied field pressure 2H Mptp and the effective pres- 
sure from the interfacial energy Act: 



Hph — 



Act 

2 Mptp' 



(32) 



where Mp and tp are the magnetization and thickness 
of the ferromagnet. When the interface is treated as 
ideally "compensated", then the exchange bias field is 
zero. On the other hand, if the AF/F interface is ideally 
uncompensated there is an interfacial energy difference 
Ac = 2Ji/a^, where Ji is the exchange coupling con- 
stant across the interface, and a is the lattice parameter 
of a simple cubic structure assigned to the AF layer. The 
EB field is ffeb = J^/a'^ Mptp (see Fig[l4| \ 

Estimating numerically the size of the EB field using 
the equation above for an ideally uncompensated inter- 
face, results in a discrepancy of several orders of mag- 
nitude with respect to the experimental observation . 
Therefore, a novel mechanism based on random fields 
at the interface acting on the AF layer is proposed as to 
drastically reduce the resulting exchange bias field. 
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FIG. 14 Schematic view of possible atomic configurations 
in a F-AF bilayer with ideal interfaces. Frustrated bonds 
are indicated by crosses. Compensated configuration a) will 
result in configuration b) by reversing the ferromagnetic spins 
through domain wall movement. It gives an exchange bias 
field of Hcb ~ Ji/ a? Mptp ■ The compensated configuration c) 
will result in the compensated configuration d). The exchange 
bias field for this case is zero {H^b ~ 0). (|84|l . 



By simple and schematic arguments Malozemoff de- 
scribes how roughness on the atomic scale of a "com- 
pensated" AF interface layer can lead to uncompensated 



^ For this example the energy is calculated as Eki = 
pair of nearest neighbor spins kl at the interface 



-JiSkSi per 



interface depicted in Fig 15 1) containing a single mono- 



atomic bump in a cubic interface gives rise to six net anti- 
ferromagnetic deviations from a perfect compensation. A 
bump shifted by one lattice spacing as shown in Fig. 15 1), 
which is equivalent to reversing the F spins, provides six 
net ferromagnetic deviations from perfect compensation. 
Thus a net energy difference of ZiJi with = 12 acts 
at the interface favoring one domain orientation over the 
other. Note that for an ideally uncompensated interface 
the energy difference is only 8Ji when reversing the F 
spins. This implies that an atomic step roughness at a 
compensated interface leads to a higher exchange bias 
field as compared to the ideally compensated interface. 

The estimates of this local field can be further refined 
assuming a more detailed model. For example, by invert- 
ing the spin in the bump shown in Fig. |15| ;, the interfacial 
energy difference is reduced by 5 x 2 at the cost of gen- 
erating one frustrated pair in the AF layer just under 
the bump. This frustrated pair increases the energy dif- 
ference by 2Ja, where Ja is the AF exchange constant. 
Thus the energy difference between the two domains be- 
comes 2 Ji + 2Ja or roughly iJ if Ji fa Ja ~ J- If one 
allows localized canting of the spins, one expects the en- 
ergy difference to be reduced somewhat further. 
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FIG. 15 Schematic side view of possible atomic moment con- 
figurations for a non-planar interface. The bump should be 
visualized on a two-dimensional interface. Configuration c) 
represents the lower energy state of a). The configuration b) 
is energetically equivalent to flipping the ferromagnetic spins 
of a). The x signs represent frustrated bonds. 



Each interface irregularity will give a local energy dif- 
ference between domains whose sign depends on the par- 
ticular location of the irregularity and whose magnitude 
is on the average 2zJ, where z is a number of order 
unity. Furthermore, for an interface which is random on 
the atomic scale, the local unidirectional interface energy 
ai — ^zJ/a^ will also be random and its average cr in a 
region will go down statistically as cr « ai/\/N, where 
N — I? I a} is the number of sites projected onto the 
interface plane. Therefore the effective AF-F exchange 
energy per unit area is given by: 



J, 



1 



eb 



N 



Jr 



1 



Ji' 



where Ji is the exchange energy of a fully uncompensated 
AF-surface. 

Given a random field provided by the interface rough- 
ness and assuming a region with a single domain of the 
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ferromagnet, it is energetically favorable for the AF to 
break up into magnetic domains, as shown schematically 
in Fig. [16] A perpendicular domain wall is the most 
preferable situation. This perpendicular domain wall is 
permanently present in the AF layer. It should be distin- 
guished from a domain wall parallel to AF/F interface, 
which according to the Mauri model develops tem- 
porarily during the rotation of the F layer. 



^^^^ 



FIG. 16 Schematic view of a vertical domain wall in the AF 
layer. It appears as an energetically favorable state of F/AF 
systems with rough interfaces H26I I84p . 



By further analyzing the stability of the magnetic do- 
mains in the presence of random fields, a characteristic 
length L of the frozen-in AF domains and their char- 
acteristic height are obtained: L « ny/ Aaf/Kaf and 
h — L/2, where Aaf is the exchange stiffness and h is 
the characteristic height of the AF domains. Once these 
domains are fixed, flipping the ferromagnetic orientation 
causes an energy change per unit area of Ac — AzJ/naL, 
which further leads to the following expression for the EB 
field: 



2z^AafKaf 



(33) 



Assuming a CoO/Co(100A) film, the calculated ex- 
change bias using the Eg. [33|is: 

^ _ 2 X 1 v/ "-",^,^-.^r-!°,= ^2.5 X 10^ [J/m^] 

7r2 X 1460 [fcA/m] 100 x 10-1° ^ 
= 580 Oe (34) 

For the estimations above we used for the exchange stiff- 
ness the following value: Aaf — Jaf/ol, where a is 

o 

the lattice parameter of CoO (a = A.27A) and Jaf = 
1.86 meV is the exchange constant for CoO (|135p . 

The characteristic length of the AF domains is for 
CoO: 



L 



tt^/Aaf/K, 



AF 



/ 0.0186 X 1.6 X 10-i9[J] 
^ ^ y 4.27 X 10-iOH X 2.5 x 10^ [J/m^] 

16.6 i (35) 



CoO(25^)/Co studied in (|132p . we notice that the cal- 
culated EB field agrees well with the value observed ex- 
perimentally. For example, the exchange bias field for 
CoO(25 i)/Co(119 A) is 557 Oe and the theoretical value 
calculated with Eq. [33] is 487 Oe. Also, the length and 
the height of the AF domains have enough space to de- 
velop. The difference between theory and experiment 
is, however, that experimentally AF domains occur after 
the very first magnetization reversal, whereas within the 
Malozemoff model the AF domains are assumed to de- 
velop during the field cooling procedure. Nevertheless, 
the agreement appears to be excellent. 



IX. DOMAIN STATE MODEL 

The Domain State model (DS) introduced by Usadel 
and coworkers (EllMllinnilinillllSlIinilllZlisa micro- 
scopic model in which disorder is introduced via magnetic 
dilution not only at the interface but also in the bulk of 
the AF layer as in Fig. [17] The key element in the model 
is that the AF layer is a diluted Ising antiferromagnet 
in an external magnetic field (DAFF) which exhibits a 
phase diagram like the one shown in Fig. 18 (I117p . In 



zero field the system undergoes a phase transition from 
a disordered, paramagnetic state to a long-range-ordered 
antiferromagnetic phase at the dilution dependent Neel 
temperature. In the low temperature region, for small 
fields, the long-range interaction phase is stable in three 
dimensions. When the field is increased at low tempera- 
ture the diluted antiferromagnet develops a domain state 
phase with a spin-glass-like behavior. The formation of 
the AF domains in the DS phase originates from the sta- 
tistical imbalance of the number of impurities of the two 
AF sublattices within any finite region of the DAFF. This 
imbalance leads to a net magnetization which couples to 
the external field. A spin reversal of the region, i. e., the 
creation of a domain, can lower the energy of the sys- 
tem. The formation of a domain wall can be minimized 
if the domain wall passes through nonmagnetic defects 
at a minimum cost of exchange energy. 



m 




FIG. 17 Sketch of the domain state model with one ferromag- 
netic layer and three diluted antiferromagnetic layers. The 
dots mark defects (|117|l . 



The height of the AF domains is h = L/2 = 8.3 A. 
Comparing this value to the experimental data on 



Nowak et al. (|117p further argue that during the field 
cooling below the irreversibility line Ti(B), in an exter- 



16 



1 



OS 

















FIG. 18 Schematic phase diagram of a three-dimentional di- 
luted antiferromagnet. AF is the antiferromagnetic phase, 
DS the domain state phase, and PM the paramagnetic 
phase ([n7|) . 



nal field and in the presence of the interfacial exchange 
field of the ferromagnet, the AF develops a frozen do- 
main state with an irreversible surplus of magnetization. 
This irreversible surplus magnetization controls then the 
exchange bias. 

The F layer is described by a classical Heisenberg 
model with the nearest-neighbor exchange constant Jp. 
The AF is modelled as a magnetically diluted Ising sys- 
tem with an easy axis parallel to that of the F. The 
Hamiltonian of the system is given by pi7p : 



magnetization tries to orient the F layer along its direc- 
tion. The coercive field has to overcome this barrier, and 
the higher the interface magnetization of the DAFF, the 
stronger is the field required to reverse the F layer. The 
interface magnetization decreases with increasing DAFF 
thickness due to a coarsening of the AF domains accom- 
panied by smoother domain walls. 




H = -Jf S,-Sj-Y.{d,Sl + d,Sl + fiBS,) 

<i,j>eF ieF 

— Jaf ^ eiej(T.i(Tj — ^ ^BztiCTi 

<i,i><LAF ieAF 

— JiNT ^ ^i<^iSjz, (36) 

<ieAF,jeF> 

where the Si and crj arc the classical spin vectors at the 
ith site of the F and AF, respectively. The first line 
contains the energy contribution of the F, the second 
line describes the diluted AF layer, and the third line 
includes the exchange coupling across the interface be- 
tween F and DAFF, where it is assumed that the Ising 
spins in the topmost layer of the DAFF interact with the 
z component of the Heisenberg spins of the F layer. 

In order to obtain the hysteresis loop of the system, the 
Hamiltonian in Eq. |36]is treated by Monte Carlo simula- 
tions. Typical hysteresis loops are shown in Fig. 19 (,101 ), 



where both the magnetization curve of the F layer and 
of the interface monolayer of the DAFF are shown. The 
coercive field extracted from the hysteresis curve depends 
on the anisotropy of the F layer, but it is also influenced 
by the DAFF. It, actually, depends on the thickness and 
anisotropy of the DAFF layer. The coercive field de- 
creases with the increasing thickness of the DAFF layer 
(jlOip . which can be understood as follows: the interface 



FIG. 19 Simulated hysteresis loops within the domain state 
model. The top hysteresis belongs to the ferromagnetic layer 
and the bottom hysteresis belongs to an AF interface mono- 
layer |PUT|> . 



The strength of the exchange bias field can be esti- 



mated from the Eq. 36 using simple ground state argu- 
ments. Assuming that all spins in the F remain parallel 
during the field reversal and some net magnetization of 
the interface layer of the DAFF remains constant during 
the reversal of the F, a simple calculation gives the usual 
estimate for the bias field (|117p : 



(37) 



where / is the number of the F layers and itiint is the 
interface magnetization of the AF per spin. iJeb is the 
notation for the exchange bias field in Ref. (|117p and is 
equivalent to H^b in this chapter. For an ideal uncom- 
pensated interface {mijfT = 1) the exchange bias is too 
high, whereas for an ideally compensated interface the 
exchange bias is zero. Within the DS model the interface 
magnetization miNT < 1 is neither a constant nor is it a 
simple quantity (I117p . Therefore, it is replaced by m/us, 
which is a measure of the irreversible domain state mag- 
netization of the DAFF interface layer and is responsible 
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for the EB field. With this, an estimate of exchange bias 
field for I = 9, Jjnt = -3.2 x 10'^^ J, and fi = 1.7 ^is 
gives a value of about 300 Oe. 

The exchange bias field depends also on the bulk prop- 
erties of the DAFF layer as shown by Miltenyi et al. (,98|l . 
There the AF layer was diluted by substituting non- 
magnetic Mg in the bulk part and away from the in- 
terface. The representative results are shown in Fig. 20 
It was shown experimentally that the EB field depends 
strongly on the dilution of the AF layer. As a function 
of concentration of the non-magnetic Mg impurities, the 
EB evolves as following: at zero dilution the exchange 
bias has finite values, whereas by increasing the Mg con- 
centration, the EB field increases first, showing a broad 
peak-like behavior, and then, when the dilution is fur- 
ther increased the EB field decreases again. Simulations 
within the DS model showed an overall good qualitative 
agreement. The peak- like behavior of the EB as a func- 
tion of the dilution is clearly seen in the simulations (see 
Fig. 20). However, it appears that at zero dilution 



the 

DS gives vanishing exchange bias whereas experimentally 
finite values are observed. The exchange bias is missing 
at low dilutions because the domains in the AF cannot be 
formed as they would cost too much energy to break the 
AF bonds. This discrepancy (|117( I122p is thought to be 
explained by other imperfections, such as grain bound- 
aries in the AF layer which is similar to dilution and 
which can also reduce the domain-wall energy, thus lead- 
ing to domain formation and EB even without dilution 
of the AF bulk. 

An important property of the kinetics of the DAFF is 
the slow relaxation of the remanent magnetization, i.e., 
the magnetization obtained after switching off the cooling 
field (117) . It is known that the remanent magnetization 
of the DS relaxes nonexponentially on extremely long- 
time scales after the field is switched off or even within 
the applied field. In the DS model the EB is related to 
this remanent magnetization. This implies a decrease of 
EB due to slow relaxation of the AF domain state. The 
reason for the training effect can be understood within 
DS model from Fig. [19] bottom panel, where it is shown 
that the hysteresis loop of the AF interface layer is not 
closed on the right hand side. This implies that the DS 
magnetization is lost partly during the hysteresis loop 
due to a rearrangement of the AF domain structure. This 
loss of magnetization clearly leads to a reduction of the 
EB. 

The blocking temperature ^ within the DS model can 
be understood by considering the phase diagram of the 
DAFF shown in Fig [TT) The frozen DS of the AF 
layer occurs after field cooling the system below the irre- 
versibility temperature Ti(h). Within this interpretation. 
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FIG. 20 In the right side of the figure is shown the fihn struc- 
ture used to study the dilution infiuence on the exchange bias 
field, a) EB field as function of the Mg concentration x in the 
Coi„a;Mga;0 layer for several temperatures, b) EB field as a 
function of difi'erent dilutions of the AF volume ()98p. 



the blocking temperature corresponds to Tiih). Since 
Ti{h) < Tjv, the blocking temperature should be always 
below the Neel temperature and should be dependent on 
the strength of the interface exchange field. The sim- 
ulations within the DS model shows that EB depends 
linearly on the temperature, as observed experimentally 
in some Co/CoO systems, but no reason is given for this 
behavior (|117p . In Ref. ([5^ the blocking temperature of 
a DAFF system (Fei_:^Zn^F2(110)/Fe/Ag with x=0.4) 
exhibits a significant enhancement with respect to the 
global ordering temperature T7v=46.9 K, of the bulk an- 
tiferromagnet Feo.6Zno.4F2. 

Overall, it is believed that strong support for the DS 
model is given by experimental observations where non- 
magnetic impurities are added to the AF layer in a sys- 
tematic and controlled fashion ([Ml Ei IMl EH 1102 HIS 
I15ip . Also, good agreement has been observed in Ref. dU, 
where the dependence of the EB as a function of AF 
thickness and temperature for IrMn/Co was analyzed. 
The asymmetry of the magnetization reversal mecha- 
nisms (jS] |6j is shown to be dependent on the angle be- 
tween the easy axis of the F and DAFF layers. It was 
found that either identical or different F reversal mecha- 
nisms (domain wall movement or coherent rotation) can 
occur as the relative orientation between the anisotropy 
axis of the F and AF is varied. This is discussed in more 
detail in sections IXILP] and IXTOI 



X. MAURI MODEL 



^ The blocking temperature of an exchange bias system is the tem- 
perature where the hysteresis loop acquires a negative or positive 
shift with respect to the field axis. It is always lower then the 
Neel temperature of the AF layer. 



The model of Mauri et al. ([53) renounces the assump- 
tion of a rigid AF layer and proposes that the AF spins 
develop a domain wall parallel to the interface. The mo- 
tivation to introduce such an hypothesis was to explore 
a possible reduction of the exchange bias field resulting 
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from the Meiklejohn and Bean model. 
The assumptions of the Mauri model are: 

• both the F and AF are in a single domain state; 

• the F layer rotates rigidly, as a whole; 

• the AF layer develops a domain wall parallel to the 
interface; 

• the AF interface layer is uncompensated (or fully 
compensated) ; 

• the AF layer has a uniaxial anisotropy; 

• the cooling field is oriented parallel to the uniaxial 
anisotropy of the AF layer; 

• the AF and F spins rotate coherently, therefore the 
Stoner-Wohlfarth model is used to describe the sys- 
tem. 




FIG. 21 Mauri model for the interface of a thin ferromagnetic 
film on a antiferromagnetic substrate (|89p . 



Schematically the spin configuration within the Mauri 
The F spins rotate coherently. 



model is shown in Fig. 21 



when the applied magnetic field is swept as to measure 
the hysteresis loop. The first interfacial AF monolayer is 
oriented away from the F spins making an angle a with 
the field cooling direction and with the anisotropy axis 
of the AF layer. The next AF monolayers are oriented 
away from the interfacial AF spins as to form a domain 
wall parallel to the interface. The spins of only one AF 
sublattice are depicted, the spins of the other sublattice 
being oppositely oriented for completing the AF order. 
At a distance ^ at the interface, a ferromagnetic layer of 
thickness tp follows. Using the Stoner-Wohlfarth model, 
the total magnetic energy can be written as (.89)) : 

E = -Ho H Mptp cos{e - P) 
+KFtF SIT? {[3) 
—Jeb cos(/3 — a) 

-2y/AAFKAF{l~COs{a)), (38) 

where the first term is the Zeeman energy of the ferro- 
magnet in an applied magnetic field, the second term is 
the anisotropy term of the F layer, the third term is the 



interfacial exchange energy and, the forth term is the en- 
ergy of the partial domain wall. The new parameter in 
the equation above is the exchange stiffness A^f- As in 
the case of the Meiklejohn and Bean model, the inter- 
facial exchange coupling parameter J^b [J/m^] is again 
undefined, assuming that it ranges between the exchange 
constant of the F layer to the exchange constant of the 
AF layer divided by an effective area (see section V.B I. 



The total magnetic energy can be written in units of 
2y/ AafKaf, which is the energy per unit surface of a 
90° domain wall in the AF layer: 



fc (1 - cos(/3)) 
-|-A [1 — COs(q; 



- fi cos(/3)^ 
/5)] + (l- 



cos(a)), (39) 



where A = J^b / {2 ^/ Aaf^af), is the interface ex- 
change, with Jeb being redefined as Jeb = ^12 /Ci 
where A12 is the interfacial exchange stiffness and 
^ is the thickness of the interface (see Fig. [21] ), 
H = Kf tp 1 1\J AafKaf is the reduced ferromagnet 
anisotropy, and k = fio H Mf tp / 2^/ AafKaf is the re- 
duced external magnetic field. 

Mauri et al. ([89|) have calculated the magnetization 
curves by numerical minimization of the reduced total 
magnetic energy Eq. 39 Several values of the A and 



H paxameters were considered providing quite realistic 
hysteresis loops. Their analysis highlights two limiting 
cases with the following expressions for the exchange bias 
field: 

-{A12/O/ ^^oMFtF /orA<l 
Heb = { (40) 
- 2^/AAFKAF/ Mo Mf tp for\->\ 

In the strong coupling limit A ^ 1 , the expression for the 
exchange bias field is similar to the value given by the 
Meiklejohn and Bean model. For this situation, practi- 
cally no important differences between the predictions 
of the two models exist. When the coupling is weak 
(A ^ 1), the Mauri model delivers a reduced exchange 
bias field which is practically independent of the inter- 
facial exchange energy. It depends on the domain wall 
energy and the parameters of the F layer. In either case 
the "1/tF" dependence is preserved by the Mauri-model. 



A. Analytical expression of exchange bias field 

In order to compare the predictions of the Mauri model 
and the Meiklejohn and Bean approach, we reconsider 
the analysis of the free energy. Starting from the ex- 
pression of the free energy Eq. |38[ the minimization with 
respect to the a and /? angle leads to the following system 
of equations: 



H 



sm 



- /3) + sin(/? - a) = 



(41) 



2 V Aaf Kaf 



sin(Q;) — sin(/3 — a) — Q 
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Similar to the Meiklejohn and Bean model we define the 



parameters P 
we set 



Kaf g^^^ jjo 



Also 



Je6 ^-^^^ ^^eb — iioMptf 

0, meaning that the applied field is swept along 
the easy axis of the AF layer. Also, we do not take into 
account the anisotropy of the ferromagnet {Kp = 0), for 
two reasons. For one the coercive fields of the exchange 
bias systems are usually much higher then the coercive 
field of the isolated F layer, and secondly it is easier to 
compare the results of the Mauri model and the M&B 
model when the anisotropy of the F layer is disregarded. 
From inspection of the equations Eq. |4l] and Eq. [23] one 
can clearly see that the first equations of the two systems 
are identical, while the second ones are different in two 
respects. The first difference is related to the term P, 
which includes the domain wall energy instead of the AF 
anisotropy term in the R ratio. The second difference 
is that instead of a sin(2 a) term in the second equation 
of Eq. 23 the Mauri model has a sin(a) term, which 
influences strongly the phase diagram shown in Fig. |22[ 
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FIG. 22 Left: The phase diagram of the exchange bias field 
and the coercive fields given by the Mauri formalism. Right: 
Typical behavior of the antiferromagnetic angle a for the two 
different regions of the phase diagram. In both regions I and 
II a shift of the hysteresis loop can exist. The coercive field 
is zero in both regimes. 



The analytical expression for the exchange bias are 
obtained by solving the system of Eqs. |4l] First step 
in solving the Eq. 41 is to extract the a angle (a = 



± arccos(± 



P-cos{3) 



=)) from the second equation 



-2 P cos{(3) ' 

and to introduce it in the first equation. Next we use the 
condition that at the coercive field (3 = — 7r/2 to obtain 
both coercive fields Hd — i?c2- Then, inserting them 
into the general expressions of Eq. [16] the coercive field 
He is zero and the exchange bias field becomes (|155p : 



H, 



eb 



Jeb '^VAaf Kaf 

Jeb P 

Mo Mf tp x/l + p2 



AF 



(42) 



This equation is plotted as a function of H/\Heb\ and 
for different P values ranging from P — Q to P = b. ( 
compare Fig. 22 left panel). The behavior of the EB field 
according to the Eq. [42] is monotonic with respect to the 
stiffness and anisotropy of the AF spins. At P ^ 1 the 
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FIG. 23 Several hysteresis loops and antiferromagnetic spin 
orientations as plotted during the magnetization reversal. For 
the simulation we used the Mauri formalism. Top row shows 
three hysteresis loops calculated for difi^erent P ratios of the 
region I shown in Fig. |22[ The right hand panel in the top 
row shows the a angle of the antiferromagnetic layer for the 
three P parameters of the hysteresis loop. The bottom raw 
are the corresponding hysteresis loops and a angles for the P 
values within region II. 



the well known expression given by M&B model. When, 
however, the P-ratio approaches low values, the exchange 
bias decreases, vanishing at P = 0, provided that the 
thickness of the AF layer is sufficiently thick to allow a 
180° wall. With some analytical analysis of the Eq. 42 
one can easily reach the limiting cases of weak coupling 
(P ^ 1) and strong coupling (P ^ 1) discussed by 
Mauri et al. ([HHj) (see Eq.jto]). In Fig. [22] right column is 



exchange bias is equal to iJ, 



which is 



shown the representative behavior of the a angle of the 
first interfacial AF monolayer as a function of the f3 ori- 
entation of the F spins during the magnetization reversal 
and for two representative P values (see the discussion 
below) . 

In Fig. 23 the hysteresis loops {m\\ = cos(/3)) and the 
corresponding AF angle rotation during the magnetiza- 
tion reversal are plotted for several P-ratios. They were 
obtained by solving numerically the system of Eqs. [41] 
For all the values of the P-ratio the magnetization curves 
are shifted to negative values of the applied magnetic 
field. We distinguish two different regions with respect 
to the behavior of the a angle of the first AF monolayer. 
In the first region , for P > 1 (region I), the AF mono- 
layer in the proximity of the F layer behaves similar to 
the Meiklejohn and Bean, namely the a angle deviates 
reversibly from the anisotropy direction as function of 
/3. The maximum value of the a angle acquired during 
the rotation of the F layer is two times higher for the 
Mauri model as compared to the M&B model, reaching 
a maximum value of 90° at P = 1. The coercive field in 
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this region is zero. The angle a in the region II where 
P < 1 has a completely different behavior. It rotates 
with the ferromagnet following the general behavior de- 
picted in Fig. |23j Notice that a follows monotonically 
the rotation of the F spins, with no jumps or hysteresis- 
like behavior in contrast to the M&B model(see Fig. |9| . 
Very importantly, the exchange bias field does not vanish 
in this region and therefore no additional coercive field 
related to the AF is observed, provided that the AF layer 
is sufficiently thick to allow for a domain wall as shown 
in Fig. 21 In this region (II) the EB field is smaller as 
compared to the M&B model. This reduction is more 
clearly seen further below, when analyzing the azimuthal 
dependence of the EB field within the Mauri model. 
Comparing the phase diagram of the Mauri model 



Fig. 22 left to the corresponding one given by the M&B 
model Fig.|9]left one can clearly see that region I of both 
models is very similar with respect to the qualitative be- 
havior of the exchange bias field as a function of the R- 
ratio and, respectively, P-ratio. However, we can com- 
pare those curves only when accounting for the variation 
of the EB field as a function of the anisotropy of the AF 
layer. Both models predict that the EB field depends 
on the anisotropy of the AF layer in a similar qualita- 
tive manner. Additionally, within the M&B model the 
EB field includes also the dependence on the thickness 
of the AF layer, which is not visible in the Mauri model. 
The other regions of both phase diagrams are completely 
different. Within the Mauri model, the exchange bias 
does not vanish at P < 1, but it continuously decreases, 
whereas the M&B model predicts that the exchange bias 
field vanishes for i? < 1 leading to enhanced coercivity. 
Also note that for the weak coupling region (II) of the 
Mauri model, the exchange bias would strongly depend 
on the temperature through the anisotropy constant of 
the AF layer (f5T]l . 



B. Azimuthal dependence of the exchange bias field 

Next we analyze the azimuthal dependence of the EB 
field by deriving an analytical expression of the EB field 
as a function of the rotation angle 9. By solving the 
second equation of the system of equation Eqs. [41] with 
respect to /?, one finds the angle a as function of (3. Using 
the condition for the coercive field as (3 — 9 — ■k/2, and 



introducing it in the first line of Eq. 41 one obtains the 
coercive fields Hd = i?c2- It follows that the coercive 
field Hc{9) = and the exchange bias field as function of 
the azimuthal angle is (|155p : 



Heb{9) = - 



J, 



eb 



2^Aaf^af cos{9) 



field is maximum parallel to the field cooling directions 
{9 = 0) only for very large P-ratios. When P approaches 
the unity, the maximum of the EB field is shifted away 
from = 0, to higher azimuthal angles and has the value: 
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eb 
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(44) 



This expression is identical to Eq. 31 of the M&B model 
and shown in Fig. [TT] The shape of the curves evolves 
from an ideal unidirectional shape at P cx3 to a skewed 
shape at P > 1. 




FIG. 24 Azimuthal dependence of exchange bias as a func- 
tion of the 6 angle. The dotted line for P = 20 can be con- 
sidered an "ideal" case. The curves are plotted according to 
Eq.pl 



In region II (Fig. 24 left) a drastic change as compared 
to region I is seen for the maximum of the exchange bias 
field as a function of azimuthal angle. Its value decreases 
monotonically towards zero according to the following 
expression: 
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eb 
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In Fig. [24] is plotted the EB bias field calculated by 
the expression above for different values of the P-ratio, 
which was also confirmed numerically. In region I the EB 



In this region the shape of the curves is also skewed for 
P-ratios close to unity, but as P decreases towards zero, 
the curves acquire a more ideal unidirectional behavior. 
Similar to region I, the maximum EB value is shifted 
away from the field cooling direction to higher azimuthal 
angles, whereas in M&B it is shifted to lower azimuthal 
angles. 

In the limit of strong R and P-ratios (P,P S> 1) the 
Mauri and M&B models give similar results. The differ- 
ences appear for the R and P-ratios which are close to 
but higher than one. This region (0 < P, P < 5) can be 
experimentally explored in order to decide in favor of one 
or the other model. In order to distinguish between the 
Mauri and the M&B model, the azimuthal dependence 
of the exchange bias offers an excellent tool because it 
is visibly different for the two models (shift to higher or 
lowe r azimuthal angles). 
sin(^)he 

presence of the planar domain wall as described by 
Mauri et al. appears to have been confirmed experimen- 
tally for first time in Ref. (IT45jl . Although the Co/NiO 
system studied by SchoU et al. did not exhibit a hys- 
teresis loop shift, the rotation angle a of the AF planar 
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wall was deduced as function of field. Also, in a recent 
publication by Gornakov et al. (l45j) show experimental 
results that are similar to the characteristic curves for 
the a angle shown in Fig. [22] An AF critical thickness 
which is often observed experimentally do not appear 
explicitly in the Mauri model. To account for the AF 
thickness dependence, Xi and White (|177p proposed a 
model which assumes a helical structure for the AF spins 
during the magnetization reversal. The temperature de- 
pendence of EB is accounted for by the Mauri model 
through the anisotropy of the AF (Reb ~ VKaf)- 
Since the anisotropy constant is rather difficult to mea- 
sure for thin films, the evidence for the predicted tem- 
perature dependence remains elusive ([3T|l . One insuffi- 
ciency of the Mauri model is the inability to predict any 
changes in coercivity. The domain wall produces only 
shifted reversible magnetization curves ([69}. Therefore, 
further refinements of the model were introduced which 
is described in the next section. 



XI. KIM-STAMPS APPROACH 
WALL 



PARTIAL DOMAIN 



The approach of Kim and Stamps dST] EH EH [701 HM 
I155P follows from the work of Neel and Mauri et al. ex- 
tending the model of an extended planar domain wall to 
the concept of a partial domain wall in the AF layer. Bi- 
quadratic (spin-flop) and bilinear coupling energies are 
used to describe exchange biased system, where the bias 
is created by the formation of partial walls in the AF 
layer. The model applies to compensated, partially com- 
pensated, and uncompensated interfaces. 

A typical hysteresis loop indicative for a partial do- 
main wall in the AF layer is shown in Fig. 25 (69). In 



saturation at /i > 0, the F spins are aligned with the ex- 
ternal field while the AF spins are in a perfect Neel state, 
coUinear with the easy axis. The interface spins are an- 
tiparallel to the F layer due to a presumably antiparallel 
coupfing. As the field is reduced and reversed, the AF 
pins the F layer by interfacial exchange coupling until 
the critical value of the reversal field is reached at he, 
where the magnetization begins to rotate. When it is en- 
ergetically more favorable to deform the AF, rather than 
breaking the interfacial coupling, a partial wall twists up 
as the F rotates. The winding and unwinding of the par- 
tial domain wall in the AF is reversible, therefore the 
magnetization is reversible (no coercivity). This mecha- 
nism is only possible if the AF is thick enough to support 
a partial wall. The magnitude of the exchange bias is 
similar to the one given by the Mauri model. Neiter the 
partial- wall theory nor the Mauri model account, how- 
ever, for the coercivity enhancement that accompanies 
the hysteresis loop shift in single domain materials, which 
is usually observed in experiments. 

The enhanced coercivity observed experimentally, is 
proposed to be related to the domain wall pinning at 
magnetic defects. The presence of an attractive domain- 
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FIG. 25 (a) Magnetization curve for the ferromag- 
net/antiferromagnet system, (b) Calculated spin structure 
at three different points of the magnetization curve. The cre- 
ation of a partial antiferromagnet domain wall can be seen in 
(iii). Only the spins close to the F/AF interface are shown. 
From Ref. ||69||. 



wall potential in the AF layer, arising from magnetic 
impurities can provide an energy barrier for domain- 
wall processes that controls coercivity. Following the 
treatment of pinning in magnetic materials by Braun 
et. al. (Uni), Kim and Stamps examined the influence 
of a pointlike impurity at an arbitrary position in the 
AF layer. As a result, the AF energy acquires, besides 
the domain wall energy, another term which depends on 
the concentration of the magnetic defects. These defects 
decrease the anisotropy locally and lead to an overall 
reduction of the AF energy. This reduction of the AF 
energy gives rise to a local energy minimum for certain 
defect positions relative to the interface. The domain 
walls can be pinned at such positions and contribute to 
the coercivity. 

Kim and Stamps argue that irreversible rotation of 
the ferromagnet due to a combination of wall pinning 
an depinning transitions, give rise to asymmetric hys- 
teresis loops. Some examples are given in Fig. 26 
The loops are calculated with an exchange defect at 
XL — 5, for three different values of defect concentration 
Pj, where denotes the defect positions in the antifer- 
romagnet, with XL = corresponding to the interface 
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FIG. 26 Defect-induced asymmetry in hysteresis loops. The 
hysteresis loops are shown for a reduced exchange defect 
at a::_L = 5 for three concentrations: (a) pj — 0.15,. (b) 
Pj — 0.45, and (c) pj — 0.75. The components of magne- 
tization parallel (M||) (dots) and perpendicular {M±) (open 
circles) to the field direction are shown. The arrows indicate 
the directions for reversal and remagnetization. The spin con- 
figuration near the interface is shown for selected field values 
below the hysteresis curves 



layer and Xl = ^af — 1 being the free surface. At low de- 
fect concentrations, the pinning potential is insufficient 
to modify partial-wall formation. The resulting magne- 
tization curve, as shown in Fig. 
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a), is reversible and 
resembles the curve obtained with the absence of impuri- 
ties. Pinning of the partial wall occurs during reversal for 
moderate concentrations, which appears as a sharp rota- 
tion of the magnetization at negative fields, as shown in 
Fig. 26 'b). During remagnetization the wall is released 
from the pinning center at a different field, thus result- 
ing in an asymmetry in the hysteresis loop. The release 
of the wall is indicated by a sharp transition in M. The 
energy barrier between wall pinning and release increases 
with defect concentration, resulting in a larger coercivity 
and reduced bias as in Fig. 26 c). 



Within this model, the asymmetry of the hysteresis 
loops is interpreted in terms of domain-wall pinning pro- 
cesses in the antiferromagnet. This explanation appears 
to be consistent with some recent work by Nikitenko et al. 
and Gornakov et al. on a NiFe/FeMn system (! 45l I106P 
who concluded that the presence of an antiferromagnetic 
wall at the interface is necessary to explain their hystere- 
sis measurements. 



XII. THE SPIN GLASS MODEL OF EXCHANGE BIAS 

To overcame the theoretical difficulties in explaining 
interconnection between the exchange bias and coerciv- 
ity, in Ref. (|133l 1 13411 is considered a magnetic state of 
the interface between F and AF layer which is magneti- 
cally disordered behaving similar to a spin glass system. 
The assumptions of the spin glass(SG) model are: 

• the F/AF interface is a frustrated spin system 
(spin-glass like) 

• frozen-in uncompensated AF spins are responsible 
for the EB shift 

• low anisotropy interfacial AF spins contribute to 
the coercivity 

Within this model, the AF layer is assumed to con- 
tain, in a first approximation, two types of AF statesv(see 
Fig. 27 1. One part has a large anisotropy with the ori- 



entation ruled by the AF spins and another part with a 
weaker anisotropy which allows some spins to rotate to- 
gether with the F spins. This interfacial part of the AF 
is a frustrated region (spin-glass-like) and gives rise to 
an increased coercivity. The presence of a low anisotropy 
AF region can be rationalized as follows: the interface 
between the F and AF layer is never perfect, therefore 
one may assume chemical intermixing, deviations from 
stoichiometry, structural inhomogeneities, low coordina- 
tion, etc, at the interface to take place. This leads to the 
formation of a transition region from the pure AF state 
to a pure F state. On average, the anisotropy of such 
an interfacial region is reduced. In addition, structural 
and magnetic roughness can provide a weak AF interface 
region. Therefore, we assume that a fraction of the frus- 
trated interfacial spins do rotate almost in phase with the 
F spins and that they mediate enhanced coercivity. We 
describe them by an effective uniaxial anisotropy Kg^J, 
because they are coupled to the presumably uniaxial AF 
layer. 

Generally, one can visualize a spin glass system (|10p as 
a collection of spins which remains in a frozen disordered 
state even at low-temperatures. In order to achieve such 
a state, two ingredients are necessary: a) there must be 
a competition among the different interactions between 
the moments, in the sense that no single configuration 
of the spins is uniquely favored by all interactions (this 
is commonly called 'frustration'); b) these interactions 
must be at least partially random. This partial random 
state will be introduced in the M&B model as an effective 
uniaxial anisotropy. 

Adding this effective anisotropy to the M&B model, 
the free energy reads ()133p : 

E = -^ioH Mptp cos{9 - (3) 
+ Kptp sin^(/3) 
-I- KAF{tAF)tAF sva^ia) 



K'/J sin2 



(/3 - 7) - J:i^ cos(/3 



(46) 
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Frozen-in AF Spins 

[Rotatable AF Spins 



FIG. 27 Schematic view of the SG Model. At the inter- 
face between the AF and the F layer the AF anisotropy is 
assumed to be reduced leading to two types of AF states af- 
ter field cooling the system; frozen-in AF spins and rotatable 
AF spins. After reversing the magnetic field, the rotatable 
spins follows the F layer rotation mediating coercivity. The 
frozen-in spins remain largely unchanged in moderate fields. 
But some of them will also deviate from the original cooling 
state. This could lead to training effects and also to an open 
loop in the right side of the hysteresis loop. At larger applied 
fields in the negative direction, the frustrated frozen-in spins 
can further reverse leading to a slowly decreasing slope of the 
hysteresis loop. A more complex antiferromagnetic state con- 
sisting of frozen magnetic domains or/and AF grains can be 
also reduced to the basic concepts depicted here. 



where, K'^q is an effective uniaxial SG anisotropy re- 
lated to the frustrated AF spins with reduced anisotropy 
at the interface, J^l^ is the reduced interfacial exchange 
energy, and 7 is the average angle of the effective SG 
anisotropy. KApitAp) is the anisotropy constant of AF 
layer. To avoid further complications for the numeri- 
cal simulations, we neglect the thickness dependence of 
the Kap anisotropy {KAp{tAF) = Kap)- We mention 
though, that this dependence could become important 
for low AF thicknesses due to finite size effects and due 
to structurally non-ideal very thin layers. From now on, 
the MCA anisotropy of the ferromagnetic layer {Kp = 0) 
will also be neglected as to highlight more clearly the 
influences of AF layer and the SG interface onto the 
general properties of the EB systems. Note that the 
Zeeman energies of the ferromagnetic-like AF interfacial 
spins are neglected in the model since they are usually 
much smaller as compared to Zeeman energy of the F 
layer. Nevertheless, they can be seen as a vertical shift 



of the hysteresis loop (frozen-in AF spins in Fig. 27 1 and 
as an additional contribution to the total magnetization 
(rotatable AF spins in Fig. 27 1. 

The model is depicted schematically in Fig. 
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the interface two rather distinct AF phases are assumed 
to occur in an EB system: the rotatable AF spins, de- 
picted as open circles and frozen-in AF spins shown as 
filled circles. After field cooling, a presumably coUinear 
arrangement is depicted in the right hand panel. After re- 
versing the magnetic field, the rotatable AF spins follow 
the F layer rotation mediating coercivity. The frozen-in 



spins remain largely unchanged in moderate fields. But 
some of them could also deviate from the original cool- 
ing state. Irreversible changes of the frustrated AF spins 
lead to training effects and also to an open loop in the 
right side of the hysteresis loop. At larger negative ap- 
plied fields, the frustrated frozen-in spins could further 
reverse leading to a slowly degreasing slope of the hystere- 
sis loop. A more complex antiferromagnetic state consist- 
ing of frozen magnetic domain state or/and AF grains 
can be also reduced to the basic two spin components 
depicted in Fig. |27j Basically, additional frozen-in spins 
can occur in the AF layer extending to the interface and, 
therefore, leading to an even more interfacial disorder. 

Next, we evaluate numerically the resulting hystere- 
sis loops and azimuthal dependence of the exchange bias 
within the SG model. When the K'qq parameter is zero, 
the system behaves ideally as described by the M&B 
model discussed in Sec lVIl the coercive field is zero and 
the exchange bias is finite. In the other case, when the in- 
terface is disordered we relate the SG effective anisotropy 
to the available interfacial coupling energy as follows: 



K^ff = (1 - ./) J, 



eb 
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eb 



(47) 



where Jeb is the total available exchange energy and / is 
a conversion factor describing the fractional order at the 
interface, with / = 1 for a perfect interface and f — 
for perfect disorder. Some basic models to calculate the 
available exchange energy were discussed in the previous 
sections. For even more complicated situations when the 
AF consists of AF grains and/or AF domains the ex- 
change energy can be further estimated as described in 
Ref. (TM). 

With these notations we write the system of equations 



resulting from the minimization of the Eq. 46 with re- 
spect to the angles a and f3: 



(1 - f) 

h sin(6l -f3)+ ^ /' sin(2 (/3 - 7)) + sin(/3 - a) = 



R sin(2 a) ~ sin(/3 - a) = 
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is the reduced applied field and 
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is the i?-ratio defining the strength of the AF layer. 

The system of equations above can easily be solved 
numerically, but it does provide simple analytical ex- 
pressions for the exchange bias. Numerical evaluation 
provides the a and (3 angles as a function of the ap- 
plied magnetic field H. The reduced longitudinal com- 
ponent of magnetization along the field axis follows from 
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TO|I = cos(P — 9) and the transverse component from 
m_L = sin(P — 9). 

With the assumptions made above the absolute value 
of the exchange bias field is directly proportional to /. 
The parameter / can be called conversion factor, as it 
describes the conversion of interfacial energy into coer- 
civity. For example, in the M&B phase diagram in the 
region II and III corresponding to reduced i?-ratios, the 
exchange bias field is zero and the coercive field is en- 
hanced as a result of such a conversion of the interfacial 
energy into coercive field. 

The idea of reduced interfacial anisotropy at the in- 
terface can be traced back to the the Neel weak ferro- 
magnetism at the surface of AF particles. Neel (|105p 
discussed the training effect as a tilting of the superficial 
magnetization of the AF domains. Later, Schlenker et 
al. (1143P suggested that successive reversals of the F mag- 
netization could lead to changes of the interface uncom- 
pensated AF magnetization and therefore provide means 
of going from one ground state to another. Such mul- 
tiple interface configurations are similar to a spin glass 
system. The spin-glass interface is further discussed by 
other authors (0 (73] [771 [571 [M] IHSl [Ml [SZl [T731). Ex- 
change bias has been observed recently for spin-glass/F 
system ([TTBJ. The interfacial magnetic disorder was ob- 
served through hysteresis loop widening below a critical 
temperature point (7). Non-coUinearity have been ob- 
served at the AF/F even in saturation (UHl [1281 1129)1 . 
The frozen spins at the interface were also observed by 
MFM (65 ) . Using element specific techniques such as soft 
x-ray resonant magnetic dichroism (XMCD) and of x-ray 
resonant magnetic scattering (XRMS), both frozen and 
rotatable AF spins can be studied (lU W[ ESI [ISl 
11311 11371 11411 1180p . The frozen-in spins appear as a shift 
of the hysteresis loop along magnetization axis, whereas 
the AF rotatable spins exhibit a hysteresis loop. More- 
over, an evidence for SG behavior is recently reported in 
thin films (|501) and AF nanoparticles (|140p . Therefore, 
we believe that there is enough experimental evidence 
to consider the interface between the AF/F layer as a 
disordered state behaving similar to a spin-glass system. 



A. Hysteresis loops as a function of the conversion factor / 

If Fig. [28] we show simulations of several hysteresis 
loops as a function of the conversion factor /. We assume 
a strong antiferromagnet in contact with a ferromagnet, 
where the interface has different degrees of disorder de- 
picted in the right column of Fig. [28] For the R ratio 
we assume the following value: R = "^'^j*^'^"^ = ^'^■^1 f 

which corresponds to a 100 A thick CoO antiferromag- 
netic layer. The field cooling direction and the measuring 
field direction are parallel to the anisotropy axis of the 
AF. The anisotropy of the ferromagnet is neglected in 
the simulations below. For the interface we have cho- 
sen a SG anisotropy oriented 10 degrees away from the 
unidirectional anisotropy orientation (7 = 10°). On the 



abscissa the reduced exchange bias field h = H/\H^\ 

{H^ = — )is plotted, which then can easily be 

compared to the M&B model. With this assumption the 



system of Eqs. 48 was solved numerically. 
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FIG. 28 Longitudinal (m|| = cos(/3)) and transverse {m±_ = 
sin(/3)) components of the magnetization for a ¥'{Kf = 
0)/AF(ii = 62.5//, 7 = lO",^ = 0) bilayer for different values 
of the conversion factor / (/ = 80 %, 60 %, 20 %). We ob- 
serve, that when the AF layer is strong [R 3> 1) the hysteresis 
loops are symmetric when measured along the field cooling 
direction. The hysteresis loops are simulated by solving nu- 
merically the Eqs. [48] In the right column is schematically 
depicted the layer structure, here the emphasis is given to 
the disorder state at the interface. The AF layer is depicted 
as consisting of magnetic domains which also contribute to 
interface disorder. 



The left column shows the longitudinal component 
of the magnetization (parallel to the measuring field 
direction) (m|| = cos(/3)) whereas the middle column 
shows the transverse component of the magnetization 
(to_l = sin(/3)). 

We observe that with decreasing conversion factor / 
the exchange bias vanishes linearly. The reduction of 
the EB field is accompanied by an increased coercivity. 
The shape of the hysteresis loop is close to the results 
found in literature. For instance the hysteresis loop with 
/ — 60% is similar to the data shown in Refs. ((44i i51)) . 
The hysteresis loop with / = 20% is similar to the data 
shown in Refs. ((71 [130). The longitudinal and transverse 
components of the magnetization show that the reversal 
mechanism is symmetric. The symmetry is directly re- 
lated to the strength of the AF layer, when no training 
effect is involved. For the examples depicted in Fig. [28] 
the R-ratio is much larger than 1 (i? ^ 1), and therefore 
the hysteresis loops are symmetric when measured along 
the field cooling direction and along anisotropy axis of 
the AF layer. The asymmetry of the hysteresis loops is 
discussed further below. 
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B. Phase diagram of exchange bias and coercive field 
within the spin glass model 

In this section we discuss the phase diagram for the 
exchange field and coercive field as a function of the R- 
value within the SG-model. The additional parameter is 
the conversion factor /. In Fig. 29 phase diagrams are 



shown for reduced exchange bias and reduced coercive 
fields as a function of the i?-ratio for four different values 
of the conversion factor /. This allows us to compare 
directly the behavior of exchange bias fields as predicted 
in the SG model and the M&B model. The reduced ex- 



change bias field plotted in Fig. 29 (left panel) is defined: 
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where the Heb is the absolute value of the exchange bias 
within the SG model and the denominator term 



^10 Mf tf 

is the exchange bias field within the ideal M&B model. 
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FIG. 29 The dependence of the reduced exchange bias field 
heb and the reduced coercive field he as a function of the R- 
ratio for four different values of the conversion factor /. 



The reduced coercive field shown in Fig. 29 is defined: 



Mo Mf tf 



where H^'^ is the absolute value of the coercive field 
within the SG model. It has no relation to the coercive 
field of the M&B model because the coercive field within 
the M&B model is considered to be a constant when the 
exchange bias is finite. 



C. Azimuthal dependence of exchange bias and coercive 
field within the spin glass model. 

In this section, the azimuthal dependence of ex- 
change bias and coercive fields within the SG model 



are discussed and compared with experimental results of 
polycrystalline /riyMnss (15 nm) /CoyoFeao (30 nm) ex- 
change bias system p33p . 

In Fig. |30] calculated magnetization components are 
plotted together with the experimental data points, and 
in Fig. [3T| 3 the azimuthal dependence of the coercive 
field and exchange bias field are plotted and compared 
to the experimental data in Fig. [31^ . The hysteresis 
loops were calculated by numerical minimizing the ex- 
pressions in Eq. |48[ The parameters used in the simula- 
tion / = 80%, R = 5.9/ f, 7 = 20° do best reproduce the 
experimental data. Furthermore, it is assumed that the 
AF layer has a uniaxial anisotropy. The MCA anisotropy 
of the F layer is neglected {Kp — 0). Therefore, the co- 
ercivity which appears in the simulations is not related 
to the F properties, but to the interfacial properties of 
the F/AF bilayer. 
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FIG. 30 Experimental (open circles)and simulated hysteresis 
loops (black lines) for different azimuthal angles. The simu- 
lated curves are calculated by the Eq. |48| with the following 
parameters: / = 80%, R = 5.9//, 7 = 20° 



First we discuss the hysteresis loops shown in Fig. [30] 
The system is cooled down in a field oriented parallel to 
the AF anisotropy direction. The hysteresis loops (solid 
lines) are simulated for different azimuthal angles 9 of the 
applied field in respect to the field cooling orientation. In 
Fig. |30| representative hysteresis curves are shown for the 
longitudinal (to||) and transverse {m±) magnetization. 
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At 9 = 0°, the magnetization curves are symmetric and 
shifted to negative fields. At = 180°, the magnetization 
curves are also symmetric but shifted to positive fields. 
At 9 = 3°, however, the longitudinal hysteresis loop be- 
comes asymmetric. The first reversal at Hd is sharp and 
the reversal at Hc2 is more rounded. This asymmetry is 
also seen in the transverse component of the magnetiza- 
tion. The F spins rotate asymmetrically: the values of 
the /3 angle depend on the external field scan direction, 
being different for swaps from negative to positive satu- 
ration as compared with swaps from positive to negative 
saturation. As the azimuthal angle increases, the coercive 
field becomes zero. For instance, a,t 9 = 20°, 90° and 160° 
there is almost no coercivity. Also, the transverse com- 
ponent of the magnetization shows that the F spins do 
not follow a 360° path, but they rotate within the 180° 
angular space. 




FIG. 31 a) The experimental coercive field (open symbols ) 
and exchange bias (filled symbols) field as a function of the 
azimuthal angle 6.d) Simulated coercive field (dotted line) 
and exchange bias field (continuous line) as a function of the 
azimuthal angle H133p . 



In Fig. 31 the coercive field and the exchange bias field 
are extracted from the experimental and simulated hys- 
teresis loops using Eq. [16] Wc distinguish the follow- 
ing characteristics of the He and H^b' the unidirectional 
behavior (« cos{9)) of the H^b as a function of the az- 
imuthal angle is (see Fig. 11 1 clearly visible; additionally, 
the behavior of the Hf,b as a function of the azimuthal an- 
gle shows sharp modulations close to the orientation of 



the AF uniaxial anisotropy; the coercive field He has a 
peak-like behavior close to the orientation of the AF uni- 
axial anisotropy, a,t 9 = 0° and 9 = 180°. In all cases we 
find an astounding agreement between calculated curves 
and experimental data. It is remarkable, that the EB 
field and the coercive field are completely reproduced by 
the SG model 

Experimentally the azimuthal dependence of the ex- 
change bias field was first explored for NiFe/CoO bilay- 
ers (3). It was suggested that the experimental results 
can be better described with a cosine series expansions, 
with odd and even terms for H^b and He , respectively, 
rather than being a simple sinusoidal function as initially 
suggested by Meiklejohn and Bean ()92j [93|l . 

The simulations shown in this section are different with 
respect to the previous reports on the angular depen- 
dence of exchange bias field (01 [39l [Ml EM EM- One 
difference is that the MCA anisotropy of the F is sup- 
posed to be negligible when compared with the coercive 
fields obtained experimentally, and the sharp features of 
the Heb are reproduced numerically rather then being de- 
scribed by cosine series expansions. Recently, Camarero 
et. al (|22p reports on very similar azimuthal dependent 
hysteresis loops as shown here. There, an elegant way 
based on asteroid curve is used to describe the intrinsic 
asymmetry of the hysteresis loops close to the 0° and 180° 
azimuthal angle. A unidirectional anisotropy displaces 
the asteroid critical curve from the origin. Therefore, 
if the applied field is not parallel to the unidirectional 
anisotropy, the field sweep line does not pass through the 
symmetry center of the asteroid critical curve leading to 
inequivalent switching fields and consequently asymmet- 
ric reversals 



D. Dependence of exchange bias field on the thicl<ness of 
the antiferromagnetic layer 

After a short inspection of the phase diagram of EB 
and coercive field (Fig. 29 1 we notice that there is a criti- 
cal value for the R-ratio at which the exchange bias van- 
ishes and the coercive field is enhanced. This critical 
value R=l depends on four parameters: the anisotropy 
of the antiferromagnet, the interfacial exchange coupling 
parameter, the thickness of the ferromagnet, and the con- 
version factor. The conversion factor further depends on 
the AF domain and/or AF grain size, if any. 

In Fig.|32|thc normalized exchange bias field is plotted 
as a function of the AF thickness and for several conver- 
sion factors. The anisotropy of the AF layer (Kaf) and 
the Jef, parameter are assumed to be constant. We no- 
tice two main characteristics of the EB field dependence 
on the AF thickness: when / has high values close to 
unity, the EB field decreases with decreasing AF thick- 
ness. However, when / is reduced, a completely different 
behavior of the EB is observed. The EB field increases as 
the thickness of the AF decreases, developing a peak-like 
feature. This peak-like behavior for the EB field at criti- 
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FIG. 32 The normalized exchange bias h^b/ f = Heb/H^ 
as a function of the tApKAp/Jeb for different values of the 
conversion factor /. The anisotropy constant Kaf and Jet. 
are assumed to be constant as the AF thickness is being varied 
for each /. An asymmetric peak like behavior of the exchange 
bias field develops for high values of the conversion factor. 



cal AF thickness is a result of enhanced coercivity which 
is accounted for by the f- factor. Also, an essential param- 
eter is the a angle, which describes the rotation of the 
AF spins during the magnetization reversal. The critical 
thickness is preserved by the SG model, but it differs in 
magnitude as compared to the M&B model. Since some 
interfacial coupling energy is dissipated as coercivity, the 
critical thickness within SG model is lower as compared 
to the corresponding one given by the M&B model. The 
critical AF thickness within the SG model follows from 
the condition R = 1: 



SG 
AF,c 



f Jeb 



r ,MB 

J ''AF,cn 



(49) 



where and t^p^^ ^^^^ ^^'^ critical thickness pre- 

dicted by the SG and M&B models, respectively. 

Experimentally the R-ratio can be tuned by chang- 
ing the thickness of the AF and keeping the other three 
parameters constant. As a result one observes a crit- 
ical thickness of the AF layer for which the EB dis- 
appears (U ni Sni [Ml [7S1 [H. This AF critical 
thickness can be qualitatively understood within the MB 
model. When the hardness of the AF layer is reduced, the 
AF spins will rotate under the torque created by the F 
layer trough the interfacial coupling constant. The shape 
of the EB as function of AF thickness , however, can be 
different from one system to another depending on the 
other three parameters. The most prominent experimen- 
tal feature of the EB dependence on the AF thickness 
is the development of a peak close to the critical thick- 
ness. Several proposals were made to describe the pecu- 
liar shapes of EB field dependence on AF thickness. Ac- 
cording to the Malozemoff model, a change in the AF do- 
main size as function of AF thickness results in a change 



of exchange bias magnitude. Another influences on the 
AF dependence of the EB and coercive field are of struc- 
tural origin [751 [75t . It has been shown by Kuch 
et. al (|78p that at the microscopic level, the coupling 
between the AF an F layers depends on the atomic layer 
filling and on the morphology of the interface. The AF-F 
coupling was observed to vary by a factor of two between 
filled and half-filled interface. Moreover, islands and va- 
cancy islands at the interface lead to a quite distinct cou- 
pling behavior. Therefore, structural configurations are 
indeed contributing to the EB-dependence as function of 
AF thickness. 

The peak-like behavior of the EB field as function of 
the AF thickness is strongly dependent on temperature. 
An almost complete set of curves, showing a monotonous 
development of the AF peak from high to low tempera- 
tures was measured by Ali et al. ([T]). The data is repro- 
duced in Fig. |33] together with the simulations based on 
DS model. Although the DS model does describe well 
some experimental observed features, some discrepancies 
still exists. For instance the development of the AF peak 
as well as the critical AF thickness as function of tem- 
perature are more pronounced in the experimental data 
as compared to the DS simulations. 
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FIG. 33 Left: IrMn thickness dependence of the exchange 
bias field He6 for a number of temperatures. Lines between 
the points are a guide to the eye. Right: Prediction of the DS 
model for the AF thickness dependence of the exchange bias 
field from the stability analysis of the interface AF domains 
at different temperatures, (from Ref. |lj)) 



The SG model, through the conversion factor f, ap- 
pears to be able to describe the evolution of the criti- 
cal AF thickness(see Fig. [32] and Fig. [33] (left)). Also 



the shape of the EB dependence on the AF thickness, 
from an almost ideal M&B type at T=290 K to a pro- 
nounced peaked curve at T=2 K is qualitatively repro- 
duced. Although not considered so far, the conversion 
factor seems to be temperature dependent. This can be 
understood if we consider the basic assumption of the SG 
model, namely the frustration at the interface. Temper- 
ature fluctuations acting on metastable spin states cause 
a variation of the SG anisotropy as function of temper- 
ature. We will not further speculate on the exact tem- 
perature dependence of the f-factor, but we mention that 
a numerical analysis for the EB dependencies shown in 
Fig. [33]would allow to untangle all the parameters in the 
R-ratio. The conversion factor is given by the shape of 
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the EB curves, fJeb can be extracted from the tempera- 
ture dependence of the EB field at high AF thicknesses, 
and finally, the anisotropy constant will be deduced from 
the critical AF thickness. Note also, that the SG model 
has the potential to even describe the different temper- 
ature dependent shapes of the EB field, namely linear 
dependence versus more rounded shape: for an AF thick- 
ness close to the critical region, the temperature depen- 
dence of the EB bias will be clearly steeper (linear-like) as 
compared to the corresponding one at higher AF thick- 
nesses (more rounded). 



E. The blocking temperature for exchange bias 

Experimentally it is found that the temperature where 
the exchange bias effect first occurs is usually lower than 
the Neel temperature of the AF layer (T/v) (109). This 
lower temperature is called blocking temperature (Is). 
For thick AF layers Tb < T/v, whereas for thin AF layers 
Tb <C T/v (|109p . Furthermore, the coercive field increases 
starting just below T/v (with some exceptions) in contrast 
to the EB field, which appears only below Tb- 

These three experimentally observed characteristic fea- 
tures can qualitatively be explained within the M&B and 
SG models. In order to have a non-vanishing EB field in 
the region with R> I, the following condition has to be 
be fulfilled: 

Kaf > 7^ = KAF,crit, (50) 
tAF 

where KAF,crit is the critical AF anisotropy for the onset 
of the EB field. For a fixed AF layer thickness, the con- 
dition above sets a critical value for the AF anisotropy 
for which the EB can exist. Considering that the AF 
anisotropy increases steadily below T/v , for large AF layer 
thicknesses the condition of Eq. [50] is fulfilled just below 
T/v, whereas for thinner AF layers this condition is ful- 
filled at a correspondingly lower temperature Tg. It is 
clear from the phase diagrams in Fig. |29] and in Fig. |32] 
that there is a region of anisotropy < Kaf < K^f 
where the EB field is zero and the coercive field is en- 
hanced. It follows that the enhancement of the coercive 
field should be observed above the blocking temperature 
and below the Neel temperature of the AF. This situa- 
tion is indeed observed experimentally. For the case of 
CoO(25 A)/Co layers the coercive field increases starting 
from the T^°^ — 291 K , whereas the exchange bias field 
first appears below Tb = 180 K (fUHll . 

Further possible causes for a reduced blocking tem- 
perature and for the behavior of the EB and coercive 
fields as a function of temperature are discussed else- 
where: finite size effects(W), stoichiometry ()168p or mul- 
tiple phases (fTM)l . AF grains ((53|) and diluted AF ([TT7)l . 



XIII. TRAINING EFFECT 

The training effect refers to the dramatic change of the 
hysteresis loop when sweeping consecutively the applied 
magnetic field of a system which is in a biased state. 
The coercive fields and the resulting exchange bias field 
versus n, where n is the n*'' measured hysteresis loop, 
displays a monotonic dependence (|1211 11421 11421 I143|l . 
The absolute value of Hd and of the EB field decreases 
from an initial value at n = 1 to an equilibrium value 
aX n — oo. The absolute value of the coercive field Hc2 , 
however, displays an opposite behavior, i.e. it increases 
with n. These features of the training effect is referred 
to as Type I by Zhang et. al. (ll82p . The other case when 
both \Hci\ and \Hc2\ decrease is called Type. II. In this 
section we deal only with the so-called Type I training 
effect. Several mechanisms were suggested as a possible 
cause of the effect. While it is widely accepted that the 
training effect is related to the unstable state of the AF 
layer and/or F/AF interface prepared by field cooling 
procedure, it is not yet well established what mechanisms 
are dominantly contributing to the training effect. 

Neel (|105p discussed the training effect as a tilting of 
the superficial magnetization of the AF domains. This 
would lead to a Type I training effect. Neel also discussed 
that a creeping effect could lead to a Type II training 
effect. 

Micromagnetic simulations within the DS model (f98l 
1117 ') show that the hysteresis curve is not closed after a 
complete loop. The lost magnetization is directly related 
to a partial loss of the superficial magnetization of the 
AF domains, which further leads to a decreased exchange 
bias. 

Zhang et. al. p8ip suggested that the training effect 
can be explained by incorporating into the Fulcomer and 
Charap's model (36) positive and negative exchange cou- 
pling between the grains constituting the AF layers. 

In Ref. ('565 ; the authors found direct evidence for the 
proportionality between the exchange bias and the to- 
tal saturation moment of the heterostructure. The find- 
ings were related to the prediction of the phenomenolog- 
ical M&B approach, where a linear dependence of the 
exchange bias on the AF interface magnetization is ex- 
pected. 

Binek (fT4| suggested that the phcnomenological ori- 
gin of the training effect is a deviation of the AF inter- 
face magnetization from its equilibrium value. Analytical 
calculations in the framework of non-equilibrium thermo- 
dynamics leads to a recursive relation accounting for the 
dependence of the H^b field on n. 

Hoffmann ([55)l argues that only biaxial AF symme- 
try can lead to training effects, reproducing important 
features of the experimental data, while simulation with 
uniaxial AF symmetry show no difference between the 
first and second hysteresis loops. 

Experiments performed by PNR, AMR and Kerr Mi- 
croscopy (Ili SZl Si Uni UMl [Sni [Ml HBJ also support 
the irreversible changes taking place at the F/AF inter- 
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FIG. 34 Simulations of training effect within the SG model. 
Longitudinal a) and transverse b) components of magnetiza- 
tion for a scenario (see text) involving irreversible changes of 
the AF domain angle 6af during the magnetization reversal. 
Longitudinal c) and transverse d) components of magnetiza- 
tion for the case when only the interface disorder parameter 
7 increases and the AF state remains unchanged. 



face and in the AF layer. It has been observed that dur- 
ing the very first reversal at Hd , interfacial magnetic do- 
mains are formed and they do not disappear even in pos- 
itive or negative " saturation" . The interfacial domains 
serve as seeds for the subsequent magnetization rever- 
sals. These ferromagnetic domains at the interface have 
to be intimately related to the AF domain state (jllip . 
Therefore, the irreversible changes of the AF domain 
state are responsible for the training effect. Furthermore 
measurements have detected out-of-plane magnetic mo- 
ments (I132j I174p hinting at the existence of perpendicular 
domain walls in the AF layer, as originally suggested by 
Malozemoff. Therefore, irreversible changes of the AF 
magnetic domains and of the interfacial domains during 
the hysteresis loops play an important role for the train- 
ing effect. 



A. Interface disorder and the training effect 

In the following we analyze the AF domains and inter- 
face contributions to the training effect. We assume that 
a gradual increase of the interfacial disorder of the F/AF 
system leads to a training effect. Within the SG model, 
the magnetic state of the F/AF interface can be mimiced 
through a unidirectional induced anisotropy K^f / , which 
is allowed to have an average direction 7, where 7 is re- 
lated to the spin disorder of the interface. Also, we will 
consider the influence of a progressive rotation of an AF 
domain anisotropy during the reversal. Both situations 



will be treated below. 

In Fig. [34| i) and b) we show first and second hysteresis 
loops (longitudinal and transverse components of mag- 
netization) calculated with the help of Eq. 48 In these 



calculations we set the conversion factor to / = 60% and 
R — 62.5/ f. We consider a drastic change of an AF 
domain which progressively rotates its anisotropy axis 
during the magnetization reversal. This situation can be 
accounted for in the SG model by replacing the a angle 
in Eq. |48]with a~ O^p, where the O^p is the orientation 
of the AF domain anisotropy. We also set the 7 angle to 
be almost zero. 

Following closely the experimental observations (|129l 
UTO,) . before the first reversal Oaf is zero, and just after 
the first reversal Oaf increases towards an equilibrium 
value. The first branch of the hysteresis loop appears 
rather sharp, therefore we assume that the AF spins 
and F spins are coUinear immediately after cooling in 
a field {Oaf = 0°). For the second branch of the 1st hys- 
teresis loop we consider that Oaf — 20°, therefore the 
second leg appears more rounded. The transition from 
Oaf = 0° to Oaf = 20° is assumed to happen right after 
or during the first reversal at Hd- This is in accordance 
with the observation that for thin CoO layers (|129l I130p 
where the disordered interface appears after the first re- 
versal at Hci- Now, the first branch of the second hys- 
teresis loop is simulated with Oaf = 20° . At the third re- 
versal, we again assume that the AF domain angle further 
increases. Therefore, the second branch of the hysteresis 
loop is simulated assuming a new value of Oaf = 30° . 

The hysteresis loops bear all the features observed ex- 
perimentally pn 1351 [T^ [T^ [Tig] [TS^ [Wl \Ym. 
More strikingly, the transverse component of magneti- 
zation shows a small step increase at Hd and larger 
increase at Hc2, on the reverse path. These are typi- 
cal features observed experimentally by PNR p29l I130|l . 
AMR dm gZl mi) and MOKE (jMl- Moreover, the trans- 
verse component at saturation behaves very close to re- 
cent observations of Brems et al. (18j) on Co/CoO bilay- 
ers and lithographically nanostructured wires. After field 
cooling and before passing through the first magnetiza- 
tion reversal in the descending field branch, the resistance 
in saturation (which is proportional to the square of the 
orientation of transverse component of magnetization) 
reaches its maximum because all spins are oriented along 
the cooling field. After going through a complete hystere- 
sis loop, the resistance at saturation is reduced, indicat- 
ing that spins in the F are rotated away from the cooling 
field. After reversing the field back to positive saturation 
the resistance does not recover its initial value. Moreover, 
the untrained state can be partially reinduced by chang- 
ing the orientation of the applied magnetic field (fT5)l 
which can be interpreted as a further indication of AF 
domain rotation during the reversal. 

Next, we consider only the interface disorder through 
a progressive change of 7 angles. In Fig. [34)3) and c) 
we show first and second hysteresis loops calculated with 
the help of Eq. 48 In these calculations we consider that 
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the AF is strong, R = 62.5/ f. For the conversion fac- 
tor we take a value of / = 60%. Also, we assume the 
average AF orientation to be parallel to the field cooling 
orientation {9af — 0). Following closely the experimen- 
tal observations, before the first reversal 7 is zero, while 
just after the first reversal, 7 increases towards an equi- 
librium value. The first branch of the hysteresis loop 
appears rather sharp, therefore we assume that the AF 
spins and F spins are coUinear immediately after cool- 
ing in a field (7 = 0). For the second branch of the 1st 
hysteresis loop we consider that 7 = 10°, therefore the 
second leg appears rather rounded. The transition from 
7 = to 7 = 10° is assumed to happen right after or 
during the first reversal at Hd- The first branch of the 
second hysteresis loop is simulated with 7 = 10°. At the 
third reversal , we again assume that the disorder of the 
interface increases. Therefore, the second branch of the 
hysteresis loop (and right after during during the reversal 
at Hci of the second loop) is simulated assuming a new 
value of 7 = 20°. 

The simulations above implies a viscosity- like behavior 
to the disordered interface (37). For example, when the 
F magnetization acquires an angle with respect to the 
unidirectional anisotropy, the torque exerted on the in- 
terfacial spins will drag them away from the initial direc- 
tion set by the field cooling. Reversing the magnetization 
back to positive directions the K^f / spins will not follow 
(completely), they remain close to this position (viscos- 
ity). This is because the maximum torque exerted by 
the F spins was already acting at negative fields, while 
for positive fields it is much reduced. When measuring 
again the hysteresis loop, at negative coercive field, the 
Keff spins will rotate even further and so on. There- 
fore, the angle of the K^ff anisotropy and/or of the AF 
domains increases after each hysteresis loop, similar to a 
rachet, causing a decreased exchange bias field. 

Comparing the hysteresis curves shown in Fig. 34 i) 



cycle index n: 



with the ones in Fig. |34p) one notices the same qualita- 
tive characteristics. This is not the case for the transverse 
magnetization curve. The F spins rotate only on the pos- 
itive side for the first case, whereas for the second case the 
F magnetization rotation covers the entire 360° angular 
range. The anisotropic magnetoresistance (AMR) and 
PNR hides the chirality of the ferromagnetic spin rota- 
tion as they provide sin^(/3) information, whereas MOKE 
is sensitive to the chirality as it provides sin(/3) informa- 
tion. Therefore, measuring both hysteresis components 
by MOKE, can help to distinguish between the domi- 
nant influence on the training effect: AF domain (or/and 
grain) rotation versus SG interface instability. The train- 
ing effect is discussed furthermore in section |XIII.B 



B. Empirical expression for the training effect 

The very first empirical expression for training 
I142p effects suggested a power law dependence of the co- 
ercive fields and the exchange bias field as a function of 
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(51) 



where k is an experimental constant. This expression 
follows well the experimental dependence of the EB field 
for n > 2, but when the very first point is included to 
the fit, then the agreement is poor. 
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FIG. 35 a) Exchange bias as function of the loop index n. 
The gray line is the best fit to the data using Eq .|51[ The open 
circles are the best fit (jl2|l to the data usi ng E g. 53 The black 

b) Temperature 
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line is the best fit to the data using Eq. 
dependence of the training effect. The lines are the best fit 
to the data using Eq. [54] 



Binek (14) has shown by using non-equilibrium ther- 
modynamics that using a recursive relation, the evolution 
of the EB field as a function of n, can be well reproduced 
for all cycle indexes (n > 1). The recursive expression 
reads: 



rrn + l 



h: 



-T{H:,^HT,Y, (52) 
where F is a physical parameter which, for n ^ 1, was di- 



rectly related to the k parameter of Eq 51 It was shown 
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that a satisfactory agreement between the Eq. 51 and 
Eq. 52 is achieved for n > 3. Therefore the approach 



of Binek appears to provide the phenomenological ori- 
gin of the hitherto unexplained power-law decay of the 
EB field with increasing loop index n > 1. The analytic 
expression (Eq. 52 1 was further tested for temperature 
dependent training effect (fr?|. 

More recently the Eq. [52] has been further refined by 
extending the free energy expansion with a correction of 
the leading term. The new equation for training effect 
reads (jTSQ]) : 



Hl\ = -7fc -i7-)^-7, iH:,-H^r, (53) 



where the new 7c parameter results from the higher order 
expansion of the free energy and hence 7c ^ 7b- The 7f, 
parameter is similar to the F in Eq. [52] Both parameters 
7b and 7c exhibits a exponential dependence on the sweep 
rate for measuring the hysteresis loop. A number of three 
fit parameters is required for both Eq. [52] and Eq. |53| 
but the last equation provides better fitting results for 
moderate sweep rates. 

In the following we analyze another type of expression 
which reproduces the dependence of the coercive field and 
exchange bias field as a function of the loop index n and 
for different temperatures. It is based on the simulations 
shown in the previous section. There it was argued that 
the training affect is related to the interfacial spin disor- 
der. With each cycle the spin disorder increases slightly, 
thereby decreasing the exchange bias field. Additional 
effects are related to the AF domain size that also affects 
the magnitude of the EB and fields. Both contri- 
butions cause a gradual decrease of exchange bias as a 
function of cycle n. They can be treated probabilisti- 
cally. We suggest the following expression to simulate 
the decrease of the EB as a function of n: 



shown: one using the empirical relation Eq. 51 the sec- 
ond one is a fit performed by Binek (il2i) using the Eq. 53 
and the third one using Eq. [54] We observe that the fit 
using the Eq. [5l] {H^ = 146.6 Oe, k = U Oe) follows 
well the experimental curve for n > 2. However, the best 
fits are obtained using Eq. [53] and Eq. [54] The best fit 
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parameters using Eq, 
7^;, = 0.00029472, and~^j, - 
eters obtained from fits to 
158 0e,A/ = 



148.632 Oe, 



-2.77210"*^. The param- 
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the data using the Eq. 
are: ii^ = lbiiUe,Af = 25.87 Oe,Pi = 4.33, = 
739.14 Oe, P, = 0.39. Within the SG approach, we distin- 
guish, indeed, a sharp contribution due to low anisotropy 
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at the interface and a much weaker decrease 
frozen-in" uncompensated spins, 
temperature dependence of the 
effect for a monocrystalline 
A)(110)/CoO(300 A)(lll)/Al203 bi- 
shown 



35 d). The 



IS siiown m Fig 
sample was field cooled in saturation to the measuring 
temperature where 31 consecutive hysteresis loops were 
measured. The fits to the data using Eq. [54] are shown 
as continues lines in Fig. 35 . 



We distinguish three main characteristics related to the 
temperature dependence of the training effect: 

• each curve shows two regimes, a fast changing one 
and a slowly decreasing tail; 

• the "relaxation times" {Pi and Pf) do not visibly 
depend on the temperature; 

• the interface transition towards the equilibrium 
state is approximatively ten times faster then the 
transition of the "frozen" spins towards their stable 
configuration. 



XIV. FURTHER CHARACTERISTICS OF EB-SYSTEMS 
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Afexp{-n/Pf) + Aiexp{-n/Pi), (54) A. Asymmetries of the hysteresis loop 



where, 77"^ is the exchange bias of the n*'* hysteresis loop, 
Af and Pf are parameters related to the change of the 
frozen spins, Ai and Pi axe parameters related to the evo- 
lution of the interfacial disorder. The A parameters have 
dimension of Oersted while the P parameters have no di- 
mension but they are similar to a relaxation time, where 
the continuous variable "time" is replaced by a discrete 
variable n. We expect that the interfacial contribution 
sharply decreases with n as the anisotropy of the inter- 
facial spins is reduced (low AF anisotropy spins), while 
the contribution from the "frozen" AF spins belonging 
to the AF domains ("frozen-in" uncompensated spins ) 
appear as a long decreasing tail as they are intimately 
embedded into a much stiffer environment. 

In the following we show fits to the " trained" exchange 
bias field. In Fig. 35 1) The EB field of thirteen con- 



secutive hysteresis loops were measured at T = 10 K 
and are plotted as a function of loop index for a 
CoO(40 A)/Fe(150 A)/Al203 bilayer. Three fits are 



A curios characteristic of EB systems is the often ob- 
served asymmetry between the two branches of the hys- 
teresis loop for descending and ascending magnetic fields. 
The hysteresis loop shape of an isolated ferromagnet is 
with no exceptions symmetric with respect to the field 
and magnetization axis. This is not the case for an ex- 
change bias system, where the unidirectional anisotropy 
and the stability of the AF can result in asymmetries 
of the hysteresis loops and of the magnetization rever- 
sal modes. One can distinguish two different classes of 
the hysteresis loop asymmetries. One of them can be 
assigned to intrinsic properties of the EB systems which 
lacks training effects, and the another one is intimately 
related to irreversible changes of the AF domain struc- 
ture during the magnetization reversal, 
i) In the first category we encounter four different situ- 
ations of asymmetryic magnetization reversal all related 
to a stable interface without training effect: 
a) the first branch of the hysteresis loop is much extended 



32 



compared to the ascending branch (see Fig. |36[ ). This 
asymmetric hysteresis loop was observed in FeNi/FeMn 
bilayers (|106p . The underlying mechanism is related to 
a Mauri type mechanism for exchange bias where a par- 
allel domain wall (exchange spring ) is formed in the AF 
layer. The reversal is understood in terms of domain wall 
pinning in the antiferromagnet (|45l 1691 1106p . 




llll 



FIG. 36 a) Asymmetric hysteresis loop (a) of a NiFe/FeMn 
bilayer and schematics of domain structure at different stages 
of magnetization reversal H1U6|I . 



b) coherent rotation during the first reversal at Hd, 
domain wall nucleation and propagation at Hc2 (see 
Fig. 37 1. Such asymmetric magnetization reversal has 
been observed by PNR in Fe/FeF2 and Fe/MnF2 sys- 
tems ([5^. This asymmetry depends on the relative ori- 



entation of the field cooling direction with respect to the 
twin structure of the AF layer. The reversal asymme- 
try mentioned above takes place when the FC is parallel 
oriented with a direction bisecting the anisotropy axes 
of the two AF structural domains. When the field cool- 
ing orientation is parallel to the anisotropy axes of one 
AF domain, the reversal mechanism is symmetric, i. e. 
for both branches of the hysteresis loop magnetization 
rotation prevails. 
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FIG. 37 a) Asymmetric hysteresis loop (b) of a Fe/MnF2 bi- 
layer and and the corresponding neutron reflectivity curves: 
a) the PNR curves recorded at the first coercive field, Hci and 
c) the PNR curve recorded at Hc2. The lack of SF reflectivity 
at Hc2 suggests that the reversal proceeds by domain wall nu- 
cleation and propagation, whereas at Hci the magnetization 
reverses by rotation (|31l I32|) . 



c) sharp reversal on the descending branch and 
rounded reversal on the ascending one (see Fig. 38 ). This 
asymmetry has recently been clarified by studies of the 
azimuthal dependence of exchange bias in IrMn/F bilay- 
ers (|22l [30l I133p . It is an intrinsic property of the EB 
bilayer systems and it takes place whenever the measur- 
ing external field is offset with respect to the field cool- 



ing orientation. By simply analyzing the geometrical as- 
teroid solutions, it becomes obvious that the sweep line 
does not symmetrically cross the shifted asteroid when 
the field is not parallel to the unidirectional anisotropy. 
Actually, this is a peculiar case of asymmetry which can 
be understood even within the phcnomcnological model 
for exchange bias (1221) which assumes a rigid AF spin 
structure. Within the SG and M&B model (133) such 
asymmetric reversals can be simulated over a wide range 
of AF thicknesses and anisotropics. Also, within the DS 
model dS] El, the effect of an offset measuring field axis 
with respect to the anisotropy axes of the AF and F layer 
can result in asymmetric reversal modes. 
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FIG. 38 Asymmetric hysteresis loop of a IrMn/CoFe bilayer 
along an offset ^ = 3° angle (|133p . 



d) the descending part is steeper, while the ascend- 
ing branch is more rounded (see Fig. 



39 1 . This asym- 



metry of the hysteresis loop needs to be distinguished 
from the previous ones, since it occurs when the external 
field is oriented parallel with respect to the anisotropy 
axis. It is observed in EB bilayers with thin antiferro- 
magnetic layers or for systems containing low anisotropy 
AF layers. We call these antiferromagnets weak antifer- 
romagnets and characterize them by the R-ratio. When 
the R-ratio is slightly higher than one (weak AF layers), 
then the asymmetry of the hysteresis loop can be repro- 
duced within the SG model. When the R-ratio is much 
higher than one (strong AF layers), then the hysteresis 
loops are symmetric as shown in Fig. [28] To account for 
the asymmetry we consider the following example where 
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FIG. 39 Longitudinal (left) and transverse components 
(right) of the magnetization vector for an week antiferromag- 
net: R= 1.1. The hysteresis loop is asymmetric: the descend- 
ing part is steeper than the ascending part. The asymmetry 
is clearly seen also in the transverse component of the mag- 
netization. 



it is essential that the AF layer is weak but the R ratio is 
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higher than 1: R = 1.1, f = 60% ,7 = 5°, and 9 = 0°. 
For these values the minimization of the free energy is 
evaluated numerically. The results are plotted in Fig. [39] 
The longitudinal and transverse components of the mag- 
netization vector is shown as a function of the reduced 
field h = H/\H^\. We clearly recognize that the hys- 
teresis loop is asymmetric: steeper on the descending leg 
and more rounded on the ascending leg. The asymmetry 
is due to the large rotation angle of the AF spins during 
the F magnetization reversal. This asymmetry has not 
received experimental recognition so far, therefore it re- 
mains a prediction of the SG model, 
ii) The second class of asymmetric hysteresis loop is di- 
rectly related to the training effect, and therefore to the 
stability of the AF layer and AF/F magnetic interface 
during the magnetization reversal. 

a)sharp reversal at Hd and more rounded reversal 
on the ascending branch, at Hc2 (see Fig. 40). Mea- 
suring subsequent hysteresis loops, the rounded rever- 
sal character does not change but appears also at H^i. 
This type of asymmetry, being related to the training 
effect, is frequently reported in different exchange bias 
systems gT] iMl [53 US] IMl Ell [Ml [Iig [1431 [Ml \M\ 
11651 11741 11811 I182p . Its underlying microscopic origin, 
has been recently demonstrated to stem in irreversible 
changes that occurs in the AF layer. Polarised Neutron 
Reflectivity measurements have revealed that at Hd the 
reversal proceeds by domain wall nucleation and prop- 
agation (gi El [HOI [Ml [Ml UMl MB ^W- At the 

second coercive field , magnetization rotation is the rever- 
sal mechanism. Moreover, by analyzing the AF/F inter- 
face (11281 1129P , it has been observed experimentally that 
a transition from a coUinear state to an non-coUinear dis- 
order state occurs. It suggests that the AF layer in (CoO 
thin layer) /F evolves from a single AF to a multiple AF 
domain state. In a AF layer that exhibits a domain state 
the anisotropy orientation in different domains is later- 
ally distributed causing a reduced coercive and exchange 
bias field. Note that a thick CoO film is suggested to be 
already in a domain state (|117p after field cooling, there- 
fore one would expect a variation of the loop asymmetry 
(and training effect) with respect to the thickness of the 
AF layer. 



B. Temperature dependence of the rotatable AF spins 
(coercivity) 
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FIG. 40 (a) MOKE hysteresis loop of a CoO/Co bilayer after 
field cooling to 50 K in an external field of 2000 Oe. The black 
dots denote the first hysteresis loop, the dotted line the second 
loop. Any further loops are not significantly different from 
the second, (b) and (c) Hysteresis loops recorded by polarized 

neutrons from the same sample but at 10 K. I , IH — 

and I — h refer to non-spin flip and spin-fiip intensities as a 
function of external magnetic field. (Ref. (|129p ') 




We discuss here the temperature dependence of the 
interfacial properties. Experimentally this can best be 
studied by an element selective method (f^ 11311 11371 
11411 1180|l to distinguish between the hysteresis of the F 
and AF layer. Element specific hysteresis loops have been 
studied for Fe/CoO (113ip . which highlights the behavior 
of the rotatable interfacial AF spins. 

The exchange bias hysteresis loops measured at the 
L3 absorption edges of Co (E=780 eV, closed symbols) 
and Fe (E=708.2 eV, open symbols) and for different 



FIG. 41 The temperature dependence of the exchange bi- 
ased hysteresis loops measured at L3 absorption edges of Co 
(E=780 eV, closed symbols) and Fe (E=708.2 eV, open sym- 
bols). Scattering angle is 261 = 32° (fT^ . 



temperatures are shown in Fig. 41 After FC to 30 K 



several hysteresis loops were measured in order to elim- 
inate training effects. Subsequentially the temperature 
was raised stepwise, from low to high T. For each tem- 
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perature an element-specific hysteresis loop at the en- 
ergies corresponding to Fe and CoO, respectively, was 
measured. The hysteresis loops of Fe as a function of 
temperature show a typical behavior. At low tempera- 
tures an increased coercive field and a shift of the hys- 
teresis loop is observed. As the temperature is increased, 
the coercive field and the exchange bias decrease until 
the blocking temperature is reached. Here, the exchange 
bias vanishes and the coercive field shows little changes 
as the temperature is further increased. 

A ferromagnetic hysteresis loop corresponding to the 
CoO layer is observed for all temperatures, following 
closely the hysteresis loop of Fe , with some notable dif- 
ferences. It appears that the ferromagnetic components 
of CoO develop higher coercive fields than Fe below the 
blocking temperature. This is an essential indication that 
the AF rotatable spins mediate coercivity between the 
AF layer and the F one, justifying the conversion factor 
introduced in the SG model. After careful analysis of 
the element specific reflectivity data (fMl I13ip . one can 
conclude that a positive exchange coupling across the 
Fe/CoO interface. The ferromagnetic moment of CoO is 
present also above the Neel temperature. Here, the AF 
layer is in a paramagnetic state, therefore the coercive 
fields for the Fe and CoO rotatable spins are equal. 



C. Vertical shift of magnetization curves (frozen AF spins) 
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FIG. 42 Element-specific Co (black) and Fe (gray) hysteresis 
loops acquired at T= 15 K after field cooling in -1-200 Oe along 
the FeF2 [001] axis, parallel to the AF spin axis. The direction 
of the cooling field and the vertical shift of the Fe loop at 
T = 15 K is indicated by arrows. Reproduced from p2Up . 



A vertical shift of magnetization has been observed fre- 
quently [5^ [s^ iMi IBS] [551 [iggi [1051 ingi [TTn ing] 

11201 11221 1140p and is considered to have several origins 
related to the different mechanisms for exchange bias. 
Within the M&B model a AF monolayer in contact to 
the F layer is assumed to be uncompensated, but still 
being part of the AF lattice. At most one could expect 
a contribution to the macroscopic or microscopic magne- 
tization equal to that of the net magnetization of an AF 
monolayer and this only by probing an AF layer consist- 
ing of an odd number of monolayers. The Mauri mecha- 



nism for exchange bias is not likely to result in a vertical 
shift of the hysteresis loop, since the AF interface is com- 
pensated. Within the SG, Malozemoff, and DS models 
for EB, a small vertical shift is intrinsic. At the inter- 
face between the AF and F layer, a number of frozen AF 
spins will be uncompensated due to the proximity of the 
F layer. Their orientation is either parallel or antipar- 
allel oriented with respect to the F spins, depending on 
the type of coupling(direct or indirect exchange). They 
contribute to the magnetization of the system. In case of 
an indirect exchange coupling mechanism, the hysteresis 
loop should be shifted downwards ([1071 1108|l . whereas in 
case of direct exchange coupling the magnetization curve 
should be shifted upwards. The other part of the inter- 
face magnetization, namely the rotatable AF spins, do 
not cause any shift since they rotate in phase with the F 
layer. Within the DS model AF domains cause an addi- 
tional shift to the macroscopic magnetization curve along 
the magnetization axis (66; 117; 122). 

To demonstrate the shift of the magnetization we dis- 
cuss a recent experiment by Ohldag ei al. ('1205 using 
XMCD (|158p . The evolution of the dichroic signal as 
function of the magnetic field, providing the element spe- 



cific hysteresis loops, is shown Fig. 42 

The sample structure is 

Pd(2 nm)/Co(2.8 nm)/FeF2(68 nm)(110)/MgF2(110)/substrate 
and has been grown by via molecular beam epitaxy. The 
F layer is a polycrystalline Co whereas the AF layer 
is a FeF2(110) untwined single crystalline layer. The 
system exhibits a positive exchange bias at large cooling 
fields ()107p . For weak cooling fields the exchange bias 
curve is shifted to the negative fields, as usual for all EB 
systems. The microscopic origin of the positive exchange 
bias is an antiferromagnetic coupling at the F/AF 
interface p07p . Although, the mechanism was clearly 
demonstrated by magnetometery measurements, the 
microscopic investigation of the AF interface provides 
more detailed information. 

The hysteresis loop of Co is typical and appears 
symmetric with respect to the magnetization axis (see 
Fig. 42). This is not the case for the AF interface mag- 
netization which shows a twofold behavior: a) some in- 
terfacial AF spins are parallel oriented with respect to 
the F spins and both are rotating almost in phase; the 
AF hysteresis loop is shifted downwards with respect to 
the magnetization axis. This seems to be a direct proof 
of the preferred antiparallel coupling between the F Co 
and the AF Fe magnetization at the AF-F interface. The 
interfacial AF Fe moments are aligned during FC by the 
exchange interaction which acts as an effective field on 
the uncompensated spins, most likely confined to the in- 
terface (|120p . An AF uncompensated magnetization in 
FeF2 can also occur due to piezomagnetism, which is al- 
lowed by symmetry in rutile-type AF compounds and 
may be induced by stresses occurring below the Neel tem- 
perature (|T2l[T3ll . 

These experimental observations can also be described 
by the SG model. During the field cooling procedure the 
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frozen-in spins depicted in Fig. 27 will be pointing op- 
posite to the F spins as to fulfill the indirect exchange 
condition at the F/AF interface. The rotatable ones re- 
main unchanged in the figure, being parallel aligned with 
the F spins. The orientation of the frozen spins is gov- 
erned by the exchange interaction at the F/AF interface, 
whereas the rotatable spins are align by the F layer and 
the external field. Upon reversal the rotatable spins will 
follow the ferromagnet whereas the frozen spins remain 
pinned, leading to a shift downwards of the interfacial 
AF hysteresis loop. As the system described above does 
not exhibits training effects, no irreversible changes oc- 
curs during the magnetization reversal. In high enough 
cooling fields, the frozen-in spins will align parallel to 
the cooling field becoming also parallel with the rotat- 
able AF and F spins. This will cause a positive shift of 
the hysteresis loop (positive exchange bias) , since the ori- 
entation of the AF frozen spins is negative with respect 
to the coupling sign (I107p . 

The experimental results described above could most 
likely be described also by DS model and the Maloze- 
moff model if magnetic domains in the AF layer will be 
confirmed. 

A relative vertical shift of the magnetization related to 
training effect is described by Hochstrat et al. (i56j) for a 
NiO(0001)/Fe(110) exchange bias system. The antifer- 
romagnet is a single crystal NiO whereas the ferromag- 
netic material is an Fe layer deposited under ultrahigh 
vacuum condition. Upon successive reversals of the F 
layer a decreased of the total magnetization was observed 
by SQUID magnetometery. The variation of the vertical 
shift as a function of the hysteresis loop index n was 
correlated to a decrease of the AF magnetization. A lin- 
ear correlation between the AF magnetization, deduced 
from the vertical shift, and the exchange bias field dur- 
ing training was found, suggesting that the training effect 
may be related to a reduced AF magnetization along the 
measuring field axis and as a function of the loop index. 



In Fig. 43 we show another situation where a vertical 
shift is visible in macroscopic magnetization curves mea- 
sured by SQUID magnetometery (|17p . The system is a 
polycrystalline CoO(2.5 nm)/Co(15 nm) bilayer grown 
by magnetron sputtering (fl^ . The hysteresis loop at 
room temperature and along the easy magnetization axis 
of the F layer is symmetric with respect to the magne- 
tization axis and shows a low coercive field. Upon field 
cooling the system through the Neel of the CoO layer 
(Tjv = 291 K) to the measuring temperature T=4.2 K 
the system is set in an exchange bias state. Then, a 
hysteresis loop is measured by sweeping the field from 
2000 Oe to -2000 Oe and back. By comparing this hys- 
teresis loop with the one measured at room temperature 
one clearly observes a vertical shift up along the mag- 
netization axis. Moreover, this up-shift is due to the 
ferromagnetic Co spins which do not fully saturate at 
-2000 Oe. Previous studies have shown that after the 
reversal at Hd, the AF CoO layer breaks into AF do- 
mains exhibiting an anisotropy distribution. Due to the 



strong coupling between the F and the AF domains, the 
F layer cannot easily be saturated. Next, after repeating 
the field cooling procedure, another hysteresis loop was 
measured by sweeping the external field to much larger 
negative value, namely to -9000 Oe. Now one observes 
that the hysteresis loop becomes more symmetric with 
respect to the magnetization axis. The saturation field, 
where the hysteresis loops closes, is about -8000 Oe. 

The example above shows that a vertical shift of mag- 
netization can be related to non a homogeneous state of 
the F layer due to non-coUinearities at the AF/F inter- 
face. The relation to the training effect is clearly seen as 
different coercive fields Hc2 depending on the strength of 
the applied fields. When a stronger field is applied in the 
negative direction, the exchange bias decreases due to a 
larger degree of irreversible changes into the AF layer. 
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FIG. 43 Hysteresis loops of a CoO(25A)/Co(180A). The 
dashed line is the hysteresis loop recorded at 300 K which 
is above the Neel temperature of the CoO layer. After field 
cooling the system to 4.2 K, a hysteresis loop is recorded 
sweeping the applied field from 9000 to -9000 Oe and back to 
9000 Oe (black line, long loop). A second loop is measured 
under the same conditions as the previous one, but the field 
range is shorter, namely between 2000 Oe to -2000 Oe and 
back to 2000 Oe. The cooling field was positive. A vertical 
shift of the hysteresis loop is clearly visible for the short loop, 
whereas the long loop appears to be centered with respect to 
magnetization axis.(Ref. (jlTj)) 



A distinct class of vertical shifts is found in exchange 
bias systems where a diluted antiferromagnet (|66l 11121 
11171 I122p acts as a pinning layer. As function of dilu- 
tion of a CoO layer by Mg impurities, as well as, partial 
oxygen pressure during the deposition, the AF layer ac- 
quires an excess magnetization seen as a vertical shift of 
the hysteresis loop. This strongly supports the fact that 
the domain state in the AF layer as well as the EB effect 
is caused and controlled the defects 1122p . 

Vertical shifts in nanostructures and nanoparticles are 
often observed due to uncompensated AF spins. For a de- 
tailed discussion we refer to recent reviews by Nogues et 
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al. (jllOp and Iglesias et al. 
on AF nanoparticles (I103( I140p 



and also recent papers 



XV. FURTHER EVIDENCE FOR SPIN-GLASS LIKE 
BEHAVIOR OBSERVED IN FINITE SIZE SYSTEMS 

Although nanoparticles are not covered in this review, 
we, nevertheless, discuss two recent instances which in 
AF nanoparticles confirm the SG behavior. One is C03O4 
nanowires p40p and the other is CoO granular struc- 
ture (|5g|l . 



A. AF Nanoparticles 
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FIG. 44 Field cooled and zero field cooled hysteresis loops 
of C03O4 nanowires. From Ref. H140|) 



Nanoparticles of antiferromagnetic materials have been 
predicted by Neel ()104p to have a small net magnetic 
moment due to an unequal number of spins on the two 
sublattices as a result of the finite size (f74|l . Hysteresis 
loops of AF nanoparticles have been observed and sev- 
eral suggestions were made to account for their weak fer- 
romagnetism (j74j) . One important finite size effect of AF 
magnetic nanoparticle is the breaking of a large number 
of exchange bonds for surface atoms. This can have a 
particularly strong effect on ionic materials, since the ex- 
change interactions are superexchange interactions. The 
deficit of exchange bonds could lead to a spin disordered 
shell exhibiting spin-glass like behavior. The last effect is 
demonstrated most recently by Salabas et al. (|140p and 
discussed further below. 

In Fig. [44|two hysteresis loops of C03O4 nanoparticles 
(8 nm diameter) prepared by a nanocasting route are 
shown (|140p . Both curves were measured at T = 2 K af- 
ter cooling the system in zero field (ZFC hysteresis loop) 
and in an external applied field of -1-4 T(FC hysteresis 
loop). Whereas the ZFC loop shows typical weak ferro- 
magnetic properties specific to AF particles, after FC a 
completely different behavior is observed. The hystere- 
sis loop is exchange biased, it is vertically shifted and 
shows training effect. Moreover, the temperature de- 
pendence (not shown) of the coercive field and exchange 
bias field increases with decreasing temperature, which 
is also a usual behavior of EB systems. More strikingly 
the FC hysteresis loop does not close on the right side 



at positive fields. This open hysteresis loop is a direct 
indication of a spin-glass like behavior, similar to the 
loops observed in pure spin glass systems (fTU|l . This be- 



havior appears to support the SG model(Fig. 27 1, where 
a reduced anisotropy is assumed to occur at the F/AF 
interface. Both frozen-in spins and rot at able spins are 
directly seen in the FC hysteresis loop. Moreover, the 
irreversible changes of the surface AF spins are causing 
the open loop. A similar open loop is also predicted by 
the DS model. Within the DS model the hysteresis loop 
of the AF interface layer is not closed on the right hand 
side because the DS magnetization is lost partly during F 
reversal due to a rearrangement of the AF domain struc- 
ture. The AF particles, however, are supposed to be sin- 
gle AF domains, therefore irreversible changes are due to 
surface effects rather than caused by AF domain kinetics. 
Certainly, at very large cooling fields, the bulk structure 
of the AF particle should be also affected. The spin-glass 
like behavior was also recently observed for cobalt ferrite 
nanoparticles (|103p . for (Mn,Fe)203_t nanograins (|123p . 
and for a eFesNCrN nanocomposite system (l38|l . 



B. Extended granular AF film 
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FIG. 45 a) ZFC hysteresis loops of [CoO [x A]/Au (60 A)] 
multilayers, b) FC hysteresis loop of a [CoO [20 A]/Au (60 A)] 
multilayer, c) ZFC and FC magnetization curves as a function 
of increasing temperature T in different fields for a [CoO [15 
A]/Au (60 A)] multilayer, d) Field dependence of Tirr raised 
to the 2/3 power for two different multilayers. (Ref. (|50p) 



Another direct experimental evidence of a spin-glass 
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like behavior was observed by Gruyters (|50p on CoO / Au 
multilayers, CoO/Cu/Fe trilayers, and CoO/Fe bilayers 
with granular structure. In Fig. 45 1) the ZFC hysteresis 
loop of a CoO/Au multilayer and for different temper- 
ature is shown. One notices that the coercive field is 
enhanced at low temperatures, but no hysteresis shift 
develops. The FC hysteresis loop (see Fig. 45 d)), how- 



ever, is almost completely shifted to one side of the field 
axis. These observations are explained by Gruyters as 
an effect of the uncompensated AF spins of the granu- 
lar film structure. The saturation magnetization of this 
granular CoO film is about 60-62 enm/gcoO which would 
result in an enormous amount of uncompensated spins 
equivalent to 22% for the observed particles volume. No 
stoichiometry influences are considered to contribute to 
this value. Although the uncompensation level is un- 
clear, the evidence of a spin-glass behavior is well demon- 
strated. In Fig. 45 :) ZFC and FC magnetization curves 
are measured as function of temperature and for differ- 
ent external fields. Two main characteristic features are 
observed for these curves: an irreversibility temperature 
Tirr, where the ZFC and FC branches of McooC^) coa- 
lesce, and a pronounced peak due to superparamagnetic 
blocking in the ZFC magnetization. The direct evidence 
of a spin-glass behavior is the field dependence of Tirr 
shown in Fig.[45}i).The existence of critical lines spanned 
by the variables temperature and magnetic field can be 
explained by mean-field theory. One of these lines has 
been predicted by de Almeida and Thouless for Ising spin 
systems (fSOIl . The Tirr raised to 2/3 power as function 
of field exhibits a linear dependence in agreement with 
the predictions Almeida and Thouless. 



C. Dependence of the exchange bias field on lateral size of 
the AF domains 

The relation between the exchange bias and the re- 
duced size effects due nano-structuring of the AF-F sys- 
tems is important from both fundamental and technolog- 
ical points of view. From a fundamental point of view, 
the reduced lateral size of both F and AF objects induces 
significant changes of the exchange bias field, coercive 
fields and also the asymmetry of the hysteresis loops dH 

I2H1 1321 ail SI 151 1571 [HI [Ml IBi [TSni in 

some systems an increased exchange bias field occurs for 
reduced lateral sizes, whereas in other cases an opposite 
behavior is reported, the exchange bias field decreases 
with decreasing the lateral length scales (|110p . We refer 
in the following to the last situation. 



In Fig. 46 are depicted schematically several lateral 
systems which are commonly used to study the influences 
of the nano-structuring onto the exchange bias proper- 
ties. A reduced lateral size of the ferromagnet (Fig. 46 j) 
and Fig. [46] ;) ) gives rise to additional shape anisotropies 
for the ferromagnet leading to a change of both coercive 
and exchange bias fields as well as a change of the hys- 
teresis shape. In Fig. 46 1) these additional anisotropies 



are minimized and therefore the dependence of exchange 
bias as function of the AF lateral size is more transpar- 
ent. For all three situations we can assume that at the 
borders defined by the geometrical nanostructures, there 
is an additional disorder extending to the interface. Even 
for the case b) where the F is nano-structured we may 
expect that the lithographic process will not always be 
stopped exactly at the interface but also affecting the AF 
layer around the dot. 

Due to nano-structuring it is natural to expect that 
at the edges of the dot there are AF spins with reduced 
anisotropy. These spins will contribute to the coercivity 
at the expense of the interfacial exchange energy. The 
effective interfacial exchange energy can be written as: 
J^^^ = FfJeb, where F is a conversion factor related to 
size effects, similar to the f defined for the interface. It 
is easy to estimate the _F-parameter, as the fraction of 
the outer shell area divided by the total dot area: F « 
A1/A2 w {n{D-d)'^)/{TTD^) = 1 -2d/D + / , where 
d is the lateral thickness of the outer shell and D is the 
diameter of the dot itself. Assuming that d « D and 
that d is constant, we obtain the general expression for 
the variation of the exchange bias field as D: 



-H. 



eb 



TjD — >oo 



1 



(55) 



where H^^~^°° is the exchange bias of an extended 
film. Within this simple model the exchange bias field 
decreases as the lateral size of the AF structure de- 
creases. This 1/D dependence is often observed experi- 
mentally (fTTUll . 



a) 





FIG. 46 A schematic view of different nanostructered sys- 
tems used to study the finite size effects on the exchange bias 
properties. 



The general behavior of the hysteresis loops as a func- 
tion of the /-parameter is depicted in Fig. |28j We see 
that for high / corresponding to a large diameter AF 
dot, the coercive field is small and the exchange bias 
is high. When, however, the lateral size of a AF ob- 
ject becomes smaller, the / ratio decreases leading to an 
increased coercive field and reduced exchange bias field. 
Transition from the top hysteresis to the bottom hystere- 
sis of Fig. [28] are usually observed due to nanostructuring 
of the AF layers (H [MlEl ITTDll . 



XVI. CONCLUSIONS 

The results presented in this chapter provides an 
overview of the physics of a F/ (Interface) /AF exchange 
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bias systems. Fundamental properties of the unidirec- 
tional anisotropy are considered and discussed. The 
Meiklejohn and Bean model as well as the Mauri model 
are considered in detail and compared, both analytically 
and numerically. The Neel approach, Malozemoff model. 
Domain State model and a model of Kim and Stamps are 
shown as they provide novel and fundamental ideas on 
the EB phenomenon. The Spin Glass model is discussed 
in even more details. Further experimental results touch- 
ing the fundamentals of exchange bias are described. 
We distinguish several outcomes of our overview: 

• The exchange bias is an interface effect, as clearly 
proved by the l/tp dependence. Deviations from 
this law were observed in the literature, but funda- 
mentally this expression is clearly well established. 

• The AF anisotropics in the bulk of the AF layer and 
at the interface to a soft ferromagnetic layer give 
rise to an impressively rich behavior of the mag- 
netic properties: the hysteresis loops can be shifted 
along both field and magnetization axis and in both 
positive and negative directions, the azimuthal de- 
pendence of exchange bias exhibits non-intuitive 
behavior such as a shift of its maximum with re- 
spect to the field cooling orientation, the hysteresis 
loops are asymmetrically shaped, etc. Some of the 
exchange bias characteristics mentioned above are 
not fully and consistently revealed experimentally 
which leaves an open window for further quests. 

• The analjrtical formulae for the exchange bias 

within M&B model depends mainly on the prop- 
erties of the AF layer, namely on its anisotropy 
and thickness. 

• The peak-like increase of the EB close to the critical 

thickness of the AF layer can be reproduced by the 
SG model. This effect is also accounted for by the 
DS model. 

• Azimuthal dependence of exchange bias could help 

to distinguish between two ideal mechanisms for 
exchange bias: Mauri versus M&B models. 

• The SG model describes the conversion of coupling 
into coercivity. 

List of acronyms 
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AF Antiferromagnet 

F Ferromagnet 

MFM Magnetic Force Microscopy 

MRAM Magnetic Random Access Memory 
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XMCD Soft X-ray Resonant Magnetic 
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XRMS Soft X-ray Resonant Magnetic 
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